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NASA's Atmospheric Effects of Aviation Project (A RAP) is developing a scientific basis for assessment of the atmospheric 
impact of subsonic and supersonic aviation. A primary goal is lo assist assessments of United Nations scientific organizations and 
hence, consideration of emissions standards by the* International Civil Aviation Organization <ICAO). Engine tests have been 
conducted at AEDC to fulfill the need of AEAfV The purpose of these tests is to obtain a comprehensive database lo be used for 
supplying critical information to the atmospheric research community. It includes ( 1 ) simulated sea-level-static test data as well 
as simulated altitude data; and (2) intrusive (extractive pro!>c) data as well as non-intrmivc (optical techniques) data. A 
ccirnmcrcMl-type bypass engine with aviation fuel was used in this test .series. The test matrix was set by parametrically selecting 
the temperature, pressure, and flow rate at sea -level -static and different altitudes to obtain a parametric set of data. 
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PREFACE 


The work reported herein was sponsored by NASA through the Atmospheric Effects of 
Aviation Pf jjcct under the Engine Exhaust Trace Chemistry Commiaee and managed through the 
NASA Lewis Research Center by Mr. Richard Niedzwiecki, NASA, and Dr. Chowen Chou Wey, 
Army Research Lab (ARL). 

The work reported herein was conducted at the Arnold Engineering Development Center 
(AEDC), Air Force Materiel Command (AFMC), Arnold AFB, TN. The work was accomplished 
by the combined efforts of Sverdrup Technology* Inc. (AEDC Group), NASA Lewis Research Cen- 
ter (NASA LeRC), Pratt & Whitney (PW), Aerodyne Research, Inc. (ARI)* and the University of 
Missouri -Rolla (UMR). Sverdrup Technology is the test contractor for AEDC, AFMC, Arnold Air 
Force Base, I'M. The AEDC efforts were performed under Job No. 2 1 00 and coordinated by the Air 
Force Project Manager, Mr. Trung Lc. In addition to being a major participant* Sverdrup provided 
technical coordination among the investigative agencies under the direction of Dr. Robert P. 
Howard* Principal Investigator and Sverdrup Project Manager. The manuscript was submitted for 
publication on July 19, P96. 

The major objective of this report was to document the results of a set of parametric measure- 
ments of engine exhaust emissions conducted on a gas turbine engine in ari AEDC simulated alti- 
tude test facility. The primary authors were: 

Dr, Robert IP, Howard 
Mr, Robert S. Hiers, Jr. 

Dr. Phillip D. Whitefield 
Dr. Donald E. Hagen 
Dr. Joda C. Wormhoudt 
Dr. Richard C. Miakc-Lye 
Mr. Richard Strange 


S verd rup/A E DC 
Sverdrup/AKDC 
University of Missouri-Rolla 
University of Missouri-Rolia 
Aerodyne Research, Inc* 

Aerodyne Research, fnc. 

Pratt and Whitney, East Hartford 


The scope of work defined in the project objectives required the services and cooperation of 
many engineers and scientists from several government and private organizations. The work and 
documentation reported herein is ari accumulation of contributions from individuals acknowledged 
below and other persons too numerous to list. Although participants were a team with contributions 
beyond individual responsibilities, primary functions associated with each include: 


NASA Project Management 


Mr. Richard Niedzwiecki NASA LeRC 
Dr. Chowen Chou Wey ARI ./NASA LeRC 
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Air Force Project Management 

Mr. Trnng V. Le 
Capt. Lee Lopes 

AF/AEDC 

AF/AEDC 

* 

Principal Investigator and 
Sverdrup Project Manager 

Dr. Robert P. Howard 

Svcrdaip/ABDC 

r 

Technical Direction 

Dr. Robert P. Howard 
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PW-Enst Hartford 
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Mr. Robert Ballard 
Mr. Noel Pujol 
Mr. Scott Grigsby 
Mr. Dave Hester 
Mr. Steve Arnold 

Sverdmp/AEDC 
Sverdrup/ AEDC 
Sverdrup/AEDC 
Sverdrup/AEDC 
Sverdrup/AEDC 


Air Force Data Release 

Mr. Tru ng V. Le 
Capt. Lee Lopes 
Lt. Roy Glassco 

AF/AEDC 

AF/AF.DC 

Wright Patterson, AFB 
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Probe Rake Design, Development, and 
Application 

Mr. Robert S. Hiers y Jr. 
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* 

NQ-UV Resonance and Continuum 
Absorption Systems 

Dr. Robert P. Howard 
Ms. Alccia Davis 
Mr. Joseph Money 

Sverdrup/AEDC 
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AR1 

Sverdrup/AEDC 


Mobile Aerosol Sampling System 

Dr. Phillip D. White field 
Dr. Donald E. Hagen 
Dr. Hnrvey Lilenfield 
Mr. Alfred R, Hopkins 
Mr. Max B. Trueblood 
Mr. Jonathan Paladino 

UMR 

UMR 

McDonnell Douglas 

UMR 

UMR 

UMR 

w 


* 


2 



AEDC'TR-96-3 


Extractive Gas Sampling 


OH Laser Absorption 


Mr. Vincent A. Zacc&rdi 
Mr. Donald G. Gardner 
Mr. Douglas J. Gonnion 
Mr. Jeffrey A. Ready 
Mr. Paul A. Jalherl 
Mr. Robert S. Hicrs, Jr. 
Dr. Robert P. Howard 
Daniel R. Catalano 

Mr. Carl W. Brasier 
Dr. James A. Drakes 
Mr. Ronald G. Porter 
Mr. Michael D. Forsythe 
Dr. Robert P. Howard 


Sverdrup/AEDC 
Sverdrup/ AEDC 
Sverdrup/AEDC 
Sverdrup/ AEDC 
Sverdrup/AEDC 
Sverdrup/AEDC 
Sverdrup/AEDC 
Sverdrup/AEDC 

Sverdrup/AEDC 

Sverdrup/AEDC 

Sverdrup/AEDC 

Sverdrup/AEDC 

Sverdrup/AEDC 


AEDC'TR-96-3 


CONTENTS 

PREFACE - 1 

1.0 INTRODUCTION 7 

2.0 TEST FACILITY AND ENGINE PARAMETERS 8 

3.0 EMISSIONS MEASUREMENT SYSTEMS 9 

3. ! Probe Rake 9 

3.2 Gas Sampling and Smoke 1^ 

3.3 Optica! Non-Intrusive Techniques 1 * 

3.4 Particulates and Aerosols 13 

4.0 SUMMARY OF RESULTS. 13 

4. i Gas Sampling Data 14 

4.2 Comparison of Optical and Gas Sampling Data 16 

4.3 Particulate and Aerosol Characterization 19 

5.0 CONCLUSIONS 20 

ILLUSTRATIONS 

fi gur e 

1 . Combustor Inlet Pressure versus Temperature for all Test Conditions 23 

2. Turbine Engine Exhaust Emissions Measurements Illustration 24 

3. Horizontal C0 3 Profile for T3 = 733 K at 9. 1 -km Altitude 25 

4. Horizontal NO Profile for T3 = 733 K at 9. 1 -km Altitude 25 

5. Horizontal EI(NO) Profile for T3 * 733 K at 9. 1 -km Altitude 26 

6. Horizontal Bar Fnel/Air Ratio Summary 26 

7. Horizontal Bar EI(NO), as N0 2t Summary for all Test Conditions 27 

8. Pressure Effect on EF(NO) 27 

9. Horizontal Bar NO* versus NO 28 

10. NO x Df/Fot, versus Engine Pressure Ratio for Late Models (1980 and Later) 28 

11. Horizontal Bar FI(CO) Summary for HI Test Conditions 29 

12. Horizontal Bar EI(THC) Summary for all Test Conditions 29 

13. Comparisons of EI(NO) at Simulated Sea-Level-Static Condition 30 

14. Comparisons of EI(NO) at 7.6-km Altitude 31 

15. Comparisons of EI(NO) at 9.1 -km Altitude 32 

36. Comparisons of Et(NO) at 12. 2-km Altitude 33 

17. Comparisons of El(NO) at 15.2-km Altitude 34 

1 8. Ratio of Emission Indices N0 2 /NO t at all Test Conditions 35 


5 


AEDC-TR-96-3 


TABLES 

Tabl e Eafie 

J . AEDC Fuel Analysis Report 36 

2. USAF Aerospace Fuels Laboratory Fuel Analysis Report 37 

3. Engine and Flight Conditions . 39 

4. Measurement Systems Used to Characterize Exhaust Emissions 41 

5. Summary' of Exhaust Emissions Characterization 42 

APPENDICES 

A. Probe Rake System 45 

B. Extractive Gas Analysis System 60 

C NO-UV Resonance And Continuum Absorption Systems 77 

D. Infrared Tunable Diode Laser System 103 

E. Aerosol Sampling Measurements (Including C0 2 Diagnostic Test) 124 

F. OH-UV Laser Absorption System 150 


6 



AED(>™*96*3 


1.0 INTRODUCTION 

The National Aeronautics Space Administration (NASA), charged with the responsibility for 
the development of advanced combustor technologies for both subsonic arid supersonic aircraft, is 
developing a database for assessment of the atmospheric impact of aircraft emissions under the 
NASA Atmospheric Effects of Aviation Project (AKAP). Quantitative levels of emissions for 
various classes of engines are vital inputs to atmospheric models developed for such assessment 
programs In support of these efforts, gaseous and paniculate emissions were measured exiting an 
engine containing an annular combustor, representative of modem combustor technology. In 
addition to emissions at altitude, data were acquired at conditions equivalent to sea-level-static to 
assess the validity of extrapolation techniques used to determine altitude emissions from measure- 
ments obtained at sea-level conditions. The combustor size was representative of a 133-1 78 IcN (30- 
40,000 Ibf) thrust commercial engine. Altitude testing was primarily at 0.8 flight Mach number for 
simulation of commercial cruise operation. The measurements were performed at Arnold 
Engineering Development Center ( AEDC) at; simulated altitude conditions ranging from sea-level- 
static (SLS) to 15.3 km. The engine was operated over a range of combustor temperatures and 
pressures at each of several altitudes to provide the atmospheric modeling community with detailed 
measureme its of altitude cruise emissions and to provide NASA a set of parametric data that will 
provide insight to exhaust emission variations as a function of combustor parameters arid altitude. 

Extractive gas and aerosol sampling and nonintrusivc optical measurement systems were 
used to characterize emissions in the exhaust. This provided independent measurement techniques 
for gaseous constituents and permitted comparison of optical nonmtmsive diagnostics with the ref- 
erence conventional gas extractive methods. Optical techniques included both ultraviolet (UV) and 
infrared (IR) for in-si tu nonintriisive measurements of gaseous constituents. Extractive techniques 
included a suite of instalments recommended by the Society of Automotive Engineers (SAE) for 
measuring gaseous species and smoke number and a unique particulate sampling system developed 
by the University of Missouri-Rolla (UM.R). The UMR system, referred to as the Mobile Aerosol 
Sampling System (MASS), is being developed in part through the NASA AHA? because the current 
industry standard for particulate characterization in jet engine emissions, the “smoke number” tech- 
nique, does not provide enough information to characterize the aerosol emitted. Although useful, 
“smoke number” does not provide the essential characteristics; total aerosol and non-volatile aero- 
sol concentration, size distribution, and hydration or growth as measured by soluble mass fraction 
and critical supersaturation spectrum. These data are essential in order to effectively assess the en- 
vironmental impact of particulate engine emissions, since they also represent the input parameters 
required by models designed to describe the production and dispersion of engine aerosol einittants 
from the combustor, through the exhaust nozzle plane ;o their eventual equilibration with the am- 
bient atmosphere. 

Exhaust samples were extracted using a multi-point cruciform rake assembly located near the 
engine nozzle exit plane with (he versatility of single-point and multiple-probe average sampling. 
Aerodynamic pressure and total temperature probes were interspersed among the gas sampling 
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probes along the horizontal arm of the cruciform rake to provide spatial measurements from which 
static pressure, static temperature, and exhaust gas velocity profiles were inferred. These measure- 
ments were taken near the centerline path of the optical measurement systems and provided flow- 
field properties used in the analyses of the optical data. 

Independent of exhaust emissions characterization, flight condition and engine performance 
parameters were monitored by test facility instrumentation. Data required for relating emissions 
measurements to test conditions and engine performance am provided in this report. 

2.0 IKS'!' FACILITY AND ENGINE PARAMETERS 

This test program measured gaseous and particulate constituents exiting an engine with a 
0.28-m : mixed flow- exhaust nozzle area and an annular combustor representative of modem 
combustors in commercial use. The nozzle was slightly convergent-divergent with an exit radius of 
30 cm. The fan (bypass) air for this engine was introduced in an annular flow surrounding the core 
How upstream of the nozzle exit without forced mixing. The engine was operated from sca-level- 
static to 15.3-krr aithuce with a range of power settings from idle to takeoff. The engine cycle is 
representative of a commercial engine in the thrust category of 133-178 kN (30-40,000 Ibf). The 
overall engine pressure ratio was approximately 30:1 . The engine was tested at a flight Mach num- 
ber near 0.8, except at the 7.6-krn altitude condition where Mach number was set near 0.7. The 
aviation fuel is specified in the accompanying fuel analysis reports prepaied by the AEDC 
Chemical Latorutory (Table 1) and the United States Air Force Aerospace Fuels Laboratory at 
Wright Patterson Air Force Base, OH (1 able 2). 

The test facility provided air directly to the engine inlet simulating total temperature and 
pressure conditions appropriate for the desired flight Mach number and altitude. This testing 
method is referred to as a direcl-conncct test configuration, llie engine exhaust vented into a 1.66- 
m-diarn diffuser connected to facility exhaust pumping stations (exhausters) that maintained 
simulated-altitude pressures within the test cell The test cell was purged with air to assist engine 
cooling, sustain desired ambient pressure, and minimize recirculation of engine exhaust gases that 
could interfere with optical measurements near the engine exit. A 7.6-cm gap between the nozzle 
exit and the diffuser inlet allowed optical access to the exhaust flow field. 

In order to protect, the engine, true sea-level-static conditions are not simulated in altitude test 
facilities. The engine inlet total pressure and temperature conditions were consistent with sea-level- 
static operation, but the ambient test cell pressure was lower than one atmosphere, as noted by a 
3.1-km altitude pressure in Table 3. For these inlet conditions, the engine performance, and thus 
exhaust emissions, should be representative of sea-level-static conditions and be fairly insensitive 
to the subatmo&pheric test cell pressure at the engine exit plane. 

The test matrix was set by parametrically varying the combustor inlet temperature (T3) or 
combustor inlet pressure (P3) at simulated SLS and four other altitudes. The T3 and P3 values 
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reported herein are actually compressor discharge stagnation temperature and pressure measure- 
ments, but closely approximate combustor inlet conditions. Facility and engine parameters for 
steady-state operating conditions are summarized in Table 3. For convenience, test conditi ns are 
noted graphically on Fig. 1 as a function of combustor pressure and temperature and grouped 
according to altitude. Facility air was preconditioned to provide ‘"standard-day” temperatures at the 
engine inlet for respective altitudes. Inlet air humidity, affected by the conditioning process, was 
monitored using a dew-point hygrometer and reported as dew-point temperature. Fuel/air ratios (by 
weight), determined from facility engine performance measurements, are given in Table 3 for both 
the hot core section and mixed hot core and fan bypass streams. Tne core fuel/air is the ratio of the 
fuel flow rate to the hot core flow rate, while the mixed fuel/air is the ratio of the fuel How rate to 
the sum of the hot core and fan bypass air flow rates. Data presented in Table 3 were acquired during 
steady-state operation o! the engine and test facility. Multiple data points acquired during steady- 
state conditions were relatively constant over time and were therefore averaged and reported as a 
single value. Clock times are included in the table to document the actual test sequence, xelate 
approximate duration per steady-state test condition, and to correlate with data from other 
measurement systems. 


3.0 EMISSIONS MEASUREMENT SYSTEMS 

A suite of optical and extractive diagnostic instrumentation was assembled for the emissions 
characterization measurements study. These are listed in Table 4 and summarized below. Details of 
the respective techniques and applications to this test program are described in the appendices. Fig- 
ure 2 illustrates relative measurement locations for the optica! techniques and probes with respect 
to the inlet of the test cell diffuser. The shaded area on the rake assembly illustrates the extent of 
the engine nozzle exit area. 

3.1 PROBE RAKE 

A. cruciform multi-point rake was used to obtain exhaust samples at individual and combined 
spatial locations allowing a mapping of exit plane core, bypass, and intermediate mixed How 
regimes. The conceptual design of the rake is illustrated schematically in Appendix A along with 
probe design details. Two distinct bars mounted 90 deg apart formed a cross as illustrated in Fig. 2. 
One was mounted vertically and one horizontally with each completely spanning the exhaust flow 
field, and tieyond. The center of the rake, and thus the horizontal bar, was 2.54 cm above the 
centerline axis of the engine nozzle. The gas sampling probes were designed to meet emissions 
certification test requirements. Since this probe rake system was to eventually be used for multiple 
tests on various engine types, the gas sampling probes were spaced equally across the rake rather 
than at positions corresponding to centers of equal area of the exhaust. The gas sampling probes on 
each radial section of cither bar could be “ganged” (valved to a common line) or sampled individ- 
ually. Individual probe sampling allowed spatial mapping of the core, bypass, intermediate partially 
mixed flow regimes, and gas recirculation regions surrounding the exhaust. “Ganged” operation 
provided spatial averages for selected regions either within or external to the exhaust How field. Ex- 
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terr.al measurements quantified concentrations of gas recirculation constituents. The vertical bar 
contained only gas sampling probes. The horizontal bar of the rake assembly contained alternating 
Mach numbcr/Flow Angularity (MFA), stagnation temperature, and gas sampling probes. Probes 
were amply spaced to minimize interferences of neighboring probe interactions with the flow field. 
The data from these probes were used to determine local static temperature, static pressure, and gas 
velocity distributions required to analyze optical data for determination of gaseous species concen- 
trations and emission indices. The probe rake assembly was designed to withstand turbine engine 
exhaust conditions over the full range of engine power settings and altitudes specified for the test 
program. The probe rake assembly was mounted inside the test cell diffuser with the longer MFA 
probes extended upstream within 7.6 cm of the nozzle exit plane, thus providing optical access for 
non-intmsive diagnostics between the nozzle exit plane and probe rake assembly. 

Two additional probes were mounted to the rake body near the center of the exhaust flow- 
field to provide exhaust gas samples to a McDonnell Douglas CO : . analyzer, discussed later, and an 
AliDC instrument not related to this measurement program. Sample lines for these probes were in- 
dependent of the gas sample distribution system, and thus provided continuous samples to these an- 
alyzers with no impact on other AEDC and UMR measurement systems. 

The test facility data acquisition systems collected data from the temperature, pressure, and 
MFA probes. These data were collected in conjunction with test facility transient and steady-state 
engine performance data. Multiple steady-state data points acquired during a steady-state engine 
condition were averaged and are reported in Appendix A for each engine condition described in 
Tabic 3. 

3.2 GAS SAMPLING AND SMOKE 

The gas sampling system is similar to systems used for commercial aircraft engine emissions 
certification testing, except for the effects of equally spaced sampling probes on spatial averages of 
the exhaust. Gas samples collected at the rake (12 cm downstream of the nozzle exit plane) were 
passed through heated lines to a sample distribution system that allowed any combination of probe 
samples to be directed to the gas analyzer systems. This allowed single- and multiple-point 
measurements within the exhaust flow field or recirculation region surrounding the exhaust flow 
field. As necessary during the test period, a stainless-steel metal bellows pump was valved in line 
between the sample distribution system and gas analyzers to raise the sample pressure slightly 
above atmospheric pressure, a requirement for gas analyzer operation. 

A suite of analyzers was used to measure concentrations of gaseous constituents which are 
reported in percent or parts per million by volume, ppmv. Mon-dispersive infrared analyzers (dry) 
were used to measure concentrations of CO and C0 2 , a chemiluminescence analyzer (dry) for 
concentrations of NO, the: same analyzer (dry) using an in-line catalytic converter for concentra- 
tions of NO*, a flame ionization detector (wet) for concentrations of total hydrocarbons (THC), and 
a paramagnetic analyzer (dry) for oxygen. Dry indicates that the samples were passed through a gas 


10 


AEOC'Tft-9&-3 


sample dryer before entering the analyzer A cold mirror dew-point hygrometer was used for real- 
time monitoring of dryer efficiency. This hygrometer was independent of the one used to monitor 
the humidity of facility air supplied to the inlet of the engine. 

The method of carbon balance was used to determine fucl/air ratio and combustion efficiency 
from the measured concentrations. This also produced estimates of exhaust H 2 C) and () 2 . Thus, 0 ; 
concentrations wore both inferred from carbon balance and measured. Also, emission indices (de- 
fined as mass of constituent per 1000 pounds of fuel) were calculated and arc presented for each 
gas sampling data point. 

An SA1£ srnoke meter was used for measurements of smoke and are ^sorted as smoke 
number, However, due to lengthy sampling timer (several minutes) require u r er smoke number 
measurement, only three were performed throughout the test period. 

Data were obtained over the full set of test conditions at altitudes ranging from equivalent 
SLS to over 15 km and engine power levels from idle to takeoff. Spatial averages were obtained 
from combinations of sampling probes within the extent of the exhaust flow field: each radial ami 
of the horizontal and vertical bars, full horizontal and vertical bars, and the full cruciform 
(horizontal and vertical bar combined). Spatial profiles were obtained by sampling individual 
probes of each bar. Real-time comparisons of vertical and horizontal species concentration profiles 
and spatial averages early in the test period confirmed symmetry and gave confidence in 
representative sampling from cither the horizontal or vertical bar. Thereafter, the horizontal bar was 
primarily assigned to the conventional gas sampling system which permitted more direct 
comparisons with data from non- intrusive optical techniques, while the vertical rake was assigned 
to the UMR-MASS (described below) for simultaneous particulate characterization. Sharing the 
rake in this manner substantially decreased data acquisition time per test condition and permitted 
data to be obtained at several test conditions that would not have been possible otherwise. 

The complete set of emissions data obtained over the test period arc presented »n Appendix 
E along with corresponding engine performance data required for data analysis. These data will be 
summarized in Section 4.0. 

3.3 OPTICAL NON-INTRUSIVE TECHNIQUES 

Optical nor* -intrusive techniques provide measurements of gaseous species concentrations 
independent of extractive probe sampling methods and thus a means for comparison of these 
techniques. These optical techniques rely on a measure of spectral transmittance (of a known wen- 
characterized radiation source) through the exhaust flow field and a characterization of the 
thermodynamic properties along the optical path. Radial profiles of static temperature, static 
pressure, and radial and axial components of flow velocities provide information for determination 
of the statistical population distribution of molecular energies line broadening and shifting, and for 
the present data, mass flux for conversion of molecular densities to emission indices. These 
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properties were inferred from probe measurements as discussed earlier Also, an analysis of optical 
transmittance data acquired at a single radial position requires an assumption on the spatial distri- 
bution of the species along the optical path, whereas radial inversion of measurements at multiple 
radial positions removes this requirement for radially symmetric flows. 

The optical measurements were performed within the 7.6-cm gap between the engine nozzle 
exit and test ceil diffuser. Figure 2 illustrates the optical beam paths relative to the entrance plane 
of the diffuser and approximate extent of the engine exhaust flow field. 

Aerodyne Research, Inc. (ARI) used coincident infrared tunable diode laser (1R-TDL) beams 
for measurements of NO, N0 2 , CO^, and H : 0 concentrations along an optical path through the 
center of the plume. The TDL system allowed simultaneous operation of two lasers with an option 
to switch one laser with a third. The interchange required only a few minutes, making data 
acquisition possible for all species daring long duration steady-state engine test conditions, but 
restricting measurements to NO, H : G, and N0 2 or CO, for many test conditions. Optical mirrors 
were used to guide the IR beam through an optical window into the test cell, along a purged path to 
the vicinity of the nozzle exit, across the exhaust flow field to a retro-reflector, and then trace a 
similar path back to the TDL system detectors. The optical window was calcium fluoride with a l- 
deg wedge to suppress interference fringes. The spectra! data, analysis process, and results are 
presented in Appendix D. 

AEDC performed spectral UV absorption measurements using both NO-resonanc? and 
continuum lamps through several, chords of the exhaust flow field and the surrounding recirculation 
gas region for determination of radial NO number density profiles. As illustrated in Fig. 2, optical 
fibers were used to conveniently transmit lamp radiation into the test cell and back out to the 
spectrometer/deieclor for each system. Inside the test cell, radiation exiting a lamp fiber was colli- 
mated, transmitted across the exhaust flow field to a mirror, and reflected back through the exhaust 
to the spectromcter/detector fiber. Line-of-sight (LOS) measurements at multiple radial positions 
were performed using slide tables to traverse the optical beam from the purged lube position to the 
centerline of the exhaust flow field. This provided spectral transmittance measurements; at radial 
positions within the exhaust and surrounding gas-recircu)ation region for each engine test 
condition. A line-by-line radiative transfer model, calibrated and verified by laboratory measure- 
ments, was used to relate measured transmittances at the second banclhcad of the NO gamma (0,0) 
band to NO number density. An “onion peel" radial inversion technique was used in conjunction 
with the radiative transfer model to determine radial NO number density profiles from the measured 
radial transmittance profiles. This provided a measure of NO density in the gas-recirculation region 
surrounding the exhaust flow field, and thus scif-corrected for molecular NO absorption along LOS 
optica! paths outside the exhaust flow field. Radial profiles of number density (cm h were convert- 
ed to profiles of fractional concentration (ppmv) and to NO emission indices. These data rind 
processes are presented in Appendix C, 
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Additionally, AEDC used a frequency-scanned UV-lnser system for measurements of OH 
number density along a fixed optical path, as illustrated in Fig. 2. The laser beam was transmitted 
into the test cell through a single-mode optica! fiber. On exiting the fiber, the beam was split with 
a portion directed to a reference detector and the remainder directed through the exhaust flow field 
to the OH absorption detector. The resolution of the laser system allowed spectra! detail of single 
OH absorption lines. These data are presented in Appendix F. 

Comprehensive data are given in respective Appendices for each optical technique. The 
emission indices for NO are summarized in Section 4.2 and compared to probe sampling data. 

3.4 PARTICULATES AND AEROSOLS 

Appendix E contains details on the particulate and aerosol characterization system em- 
ployed in this study-thc University of Missouri-Rolla Mobile Aerosol Sampling System 
(MASS). A brief overview of the system is given here. For this test, the MASS was divided into 
two work stations. The first, located close to the extractive sampling probe selection manifolds, 
was used in conjunction with the extractive sampling system to acquire total condensation nuclei 
(CN) concentrations for all aerosols in real time and to fill and pressurize sample tanks. Exhaust 
samples were delivered to the MASS from individual and ganged (multiple) probe sources. An 
isobaric cooler for supersaturation pulse simulation, a laser aerosol spectrometer tor real-time 
large diameter aerosol measurements (diameters from 1 .0 to 30 (pm)), and a needle-to -grid elec- 
tric precipitator for aerosol collection on electron microscope grids, were also located at the first 
work station. The second work station consisted of two Electrostatic Aerosol Classifiers (EAC's) 
configured in series. With this system the tank samples taken at the first work station were ana- 
lyzed for size distribution and growth and/or hydration (critical supernatural ion) information for 
aerosols in the submicron range, 1 0 Tim to 500 nm. Comprehensive data presented in Appendix 
E are summarized in Section 4.3. 

4.0 SUMMARY OF RESULTS 

As mentioned earlier, the engine bad a mixed flow exhaust with the hot core and fan bypass 
air exiting the same tailpipe. The riegre** of mixing of the two flow streams was not known before 
the test. Tabic 5 summarizes the data obtained from each of the conventional extractive gas 
analyzers at each engine set point. The gas sampling data presented in the summary table are spatial 
averages measured on the horizontal bar of the rake system, which more directly compare to optical 
measurements. The summarized MASS data are typical of particle and aerosol emission measure- 
ments for each engine condition and were obtained primarily on the vertical bar of the rake assem- 
bly. This section of the report provides an overview of the results. The reader, again, is directed to 
respective appendices for more complete data and detailed discussions. 
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4*1 GAS SAMPLING DATA 

Figure 3 shows the spatial C0 2 concentration (percent by volume) distribution measured at 
individual probes along the horizontal bar for T3 ~ 733 K at the 9. 1 -km altitude, a condition typical 
for a r engine on a commercial aircraft operating at cruise conditions. Radial positions arc normal- 
ized to the nozzle exit radius, Re = 30 cm, so that the nozzle centerline is located at R/Re - 0. The 
fuel burned on a molar basis is the sum of the moles of the species containing all of the carbon 
atoms. Thus, neglecting carbon emitted in the form of particulates and trace chemical species, the 
CO; concentration is directly related to the local fuel/air ratio. As indicated in Fig. 3, C0 2 
concentrations are high over the central region of the exhaust, comprised mainly of hot core flow, 
and drop off toward the outside where the fan air mixes with the core flow. At the outer edges of 
the exhaust, the flow is predominantly cooler fan bypass air. This is consistent with what one 
expects fora mixed flow' exhaust with no forced mixing. The slight decrease in C0 2 at the center of 
the flow is expected based on the thermal design profile of the combustor exit. The radial profile 
data from the horizontal bar show good spatial symmetry about the centerline of the nozzle. 

Figure 4 is the corresponding NO concentration profile for the same engine condition. The 
shape* is similar to the CO, profile and again exhibiting good radial symmetry. The NO emission 
index, EI(NO), is proportional to the ratio of the NO concentration and fuel flow rate (C0 2 
concentration). The radial profile of EI(NO), reported as NG 2v is presented in Fig. 5. The FJ(NO) 
profile is flat across the entire radial extent of the exhaust. This behavior is consistent with a 
uniform engine core flow'. Most of the variation in both NO and C0 2 concentration profiles is 
toward the outermost regions of the tailpipe radius, and this variation is attributed to dilution of the 
core flow with fan bypass air. The limited scatter exhibited in the data indicates that both the engine 
and gas sampling system were stable and did not exhibit temporal variations for the duration of an 
engine test condition period. This stability was further substantiated by small variations in repeated 
data points during selected engine test conditions throughout the test program. 

The gas sampling profile data obtained on the horizontal and vertical bars for this and other 
engine test conditions arc presented in Appendix B. Profiles for each constituent measured by the 
gas sampling system were radially symmetric and well behaved. This indicates that the gas 
sampling data could serve as a valid reference for comparisons of the species concentrations 
determined from optical techniques. 

Puel/air ratios (by weight) were calculated for each gas sample data point using the method 
of carbon balance. The fuel carben-to-hydrogen ratio used in the calculations was obtained from 
the fuel analysis. The sampling (carbon balance) fuel/air ratios arc given in Table 5 for horizontal 
bar averages corresponding to each test condition, and arc compared to the engine performance core 
and mixed fnel/air ratios in Fig. 6. As expected, the sample fuel/air ratios arc lower than the core 
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fuel/air and greater than mixed fuel/air ratios. A sampling bar with equally spaced probes Prather 
than at positions of equal area) across the extent of the exhaust of a mixed flow engine should give 
a fxiel/air ratio lower than the core fucl/air ratio due to sampling of some fan bypass air. Also, the 
sampling fucl/air ratio should be higher than the fully mixed fuel/air ratio since the equally spaced 
samples bias the average toward the core. The sampling fuel/air ratios do fall closer to the core fuel/ 
air than the fully mixed f'uel/air. Figure 6 comparisons again indicate that the engine and 
instrumentation operated in a consistent manner throughout the test period. 

Figure 7 summarizes El(NO) versus combustor inlet temperature (T3) for all test conditions, 
again presenting only horizontal bar averages. Of interest is the variation of EI(NO) as a function of 
altitude. At each altitude, die BI(NO) increases monotonically with increasing T3, while for any 
given T3, El(NO) decreases with increasing altitude (decreasing P3). Figure 8 shows the pressure 
effect on EI(NO) for fines of constant 73. Preliminary curve fits of these data indicate that atconstant 
T3, EI(KO) varies as P3 raised to a power in the range of 0.3 to 0.5. These trends are consistent with 
existing correlations based on theoretical and rig data measured on other test programs. 

Figure 9 is a correlation of NO, versus NO concentrations for all engine test conditions. 
Variation iri the concentration ratio NO/NG t is small. For all data measured. NO concentrations 
were on average 93 percent of the NO x . It should be noted here that while both NO and NO, were 
measured by the same instalment, measurements were always sequential and the sample passed 
through a catalytic converter for measurements of NO x . The small amount of scatter about the linear 
curve fit of the data shown in Fig. 9 is indicative of the high level of engine and instrumentation 
stability at each engine test condition. 

The International Civil Aviation Organization (ICAO) requires, for purposes of regulatory 
compliance of commercial engines, the measurement and reporting of gaseous pollutants in terms 
of Q JV /F 0 „. The proper calculation and parameters involved are described in Volume II of the 
International Standards and Recommended Practices, “Environmental Protection/’ Annex 16 to 
the Convention of International Civil Aviation. D r is the mass of the gaseous pollutant emitted 
during the reference emissions landing and take-off (LTO) cycle. F {)0 is the maximum power/ 
thrust available for take-off under normal operating conditions at International Standard Atmo- 
sphere (as defined in the manual, Document 7488) SLS conditions. These are the atmospheric 
conditions to which all engine performance data should be corrected. The sco-level-static NO, 
data were corrected to Annex 16 reference specifications and applied to the performance condi- 
tions and LTO cycle of a. midsize commercial engine of the 133-kN (30,000 Ibf) thrust class. The 
subsequently calculated NO, for a one engine characteristic value (g/)<N) is shown in Fig. 
10 in relation io other late model engines. This calculation was for a pseudo-comparison of this 
engine with other engines operating in the current fleet. Data for a true were beyond the 

scope of this test program. 
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Figure 1 1 shows Ef(CO) determined from horizontal bar measurements of CO concentrations 
for all engine test conditions. The EI(CO> levels are highest for idle and off-idle engine power 
settings, represented by the lowest T3 values in the figure. Except for the idle conditions, El(CO) 
levels were very low. For fixed T3 settings, EI(CO) appears to increase with altitude which is 
especially evident from data at the highest altitude (15.2 km). 

Figure 12 summarizes the EI(THC) for all engine test conditions. Tn general the data are 
within the zero scatter (range of readings for clean air) of the instrument. At the lowest T3 there is 
an indication that El(THC) increased slightly, but the concentration levels were too low to make 
definitive statements. 

The smoke meter, recommended by ‘•'AE, was used to measure smoke from the engine during 
three high power points. The data arc presented in the summary sheet. Table 5. Smoke was basically 
not measurable. All smoke numbers were less than 1 .5, which is within the zero scatter (range of 
readings for clean air) of this technique. 

4.2 COMPARISON OF OPTICAL AND GAS SAMPLING DATA 

Species number densities, and subsequently concentrations, were deduced from the optical 
technique measurements. Thermodynamic properties along the optical paths were provided by 
probe measurements as discussed earlier. Further, emission indices we re determined for respective 
species concentrations using exhaust velocity profiles inferred from probe measurements and fuel 
flow rates from test facility engine performance data. 

UV resonance arid continuum absorption measurements were performed at all engine test 
conditions. The resonance absorption system provided a measure of NO density over the full range 
of test conditions, whereas the less sensitive continuum system did not exhibit absorption features 
at test conditions with lowest NO concentrations. Preliminary analysis of the continuum absorption 
data indicated genera! agreement with resonance absorption results, but greater uncertainty. A 
complete analysis of the continuum absorption data will be published at a later date. The following 
discussions will be restricted to UV resonance absorption data. 

UV resonance absorption measurements at multiple lines of sight allowed NO density 
determination without assuming spatial NO concentration profiles l>efore (he analysis. Radial inver- 
sions performed on the optical UV measurements through several chords of the exhaust flow and sur- 
rounding recirculation gas gave radial NO number density (cm- , ) profiles for each engine test 
condition, For more direct c omparisons with the gas sampling data profiles, these densities were con- 
verted to volumetric fractions (pptnv) us ; ng static flow-field properties inferred from probe rake data 
and the ideal gas law. Comparisons with probe sampling data (Appendix C) show that NO concentra- 
tion profiles determined by UV measurements within the exhaust region are similar in shape but with 
somewhat lower peak levels. Greatest differences were observed at simulated sea-level-static engine 
test conditions. UV and probe sampling NO concentrations agreed well for the gas-rccirculation 
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region, both showing appreciably greater NO concentrations (relative to exhaust levels) for the sea- 
level- static condition. 

NCMJV concentration profile comparisons were limited to the few test conditions for which 
probe sampling radial profile data were available. For comparisons over all engine test conditions, 
emission indices were calculated from NO-UV resonance absorption concentration profiles using 
exhaust flow' velocity profiles derived from the rake data and fuel flow rates measured by test facil- 
ity instrumentation. El(NO) data, reported as equivalent NO* are presented in Table 5. UV E:I(NO) 
data are in good agreement with probe sampling data, matching well within respective uncertainty 
bands. However, NO-UV EI(NO) tended to be lower (but not consistently) throughout the test pe- 
riod. 


An engineering statistical analysis process was applied to determine confidence bands with 
respect to measurement precision. The variance-covariance matrix developed for a set of radial 
profile transmittance*. took into account codependency on lamp reference spectra, curve-fit 
correction factors, and the linear smoothing algorithm. This process produced one-sigma variances 
for NO number density at each radial location per set of data. Overall, uncertainly in NO number 
density per radial location for several test conditions investigated ranged from ±6 to ±18 percent, 
with the greater uncertainties observed for data acquired during higher engine power settings at the 
simulated SLS condition. Statistical uncertainty at radial positions in the surrounding gas- 
recirculation region ranged from ±15 to ±25 percent Extending the statistical analysis to emission 
indices calculations, uncertainty ranged from ±7 to ±13 percent when not accounting for 
uncertainty in exhaust radial velocity profiles. It should be noted that absorption in the (0,0) band 
was fairly insensitive to static temperature and pressure, and thus contributed little to the 
uncertainty in NO number density. However, uncertainty in velocity, especially near the outer edge 
where abrupt changes; were not quantifiable from the sparse rake measurements, increased the 
reported uncertainty in emission indices. 

This test program has provided the first comprehensive set of simultaneous UV resonance 
absorption and probe sampling data for direct comparisons of NO concentrations since extensive 
studies performed in the late 1970's, The results reported here substantiate for a wide range of 
engine exhaust conditions that agreement will be achieved if probes are designed properly for the 
exhaust flow-field conditions, and proper measurement procedures are followed. Additionally, UV 
resonance absorption measurements performed in this test program demonstrate the maturity and 
robustness of this measurement system for turbine exhaust applications in harsh environments of 
altitude engine test facilities, Also, the use of fiber optics, a relatively simple traversing system, and 
linear array detectors significantly decreases the system complexity and increases system reliability 
for measurement applications. 

The NO emission index values, given in Table S for the three techniques, 1R, UV, and gas 
sampling, arc plotted in Figs. 13 through 17 for respective altitudes and displayed as functions of 
combustor inlet temperature, T3. (In these figures, the dotted lines simply serve to connect points 
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and indicate changes in ordering of the three values.) The considerations involved in assigning 
uncertainty bars for each diagnostic technique are discussed in the appendices devoted to each 
measurement method. The widths of the uncertainty bar caps on the plots are different to help 
distinguish which uncertainty applies to each technique: the width is least for the gas sampling 
EI(NO), and largest for the IR diode laser absorption EI(NO). 

All three techniques display similar trends and agree within their uncertainty limits. No 
systematic ordering among (he values from the three methods is seen, even at levels within the error 
limits. In Fig. 1 4, BI(NO) for the higher engine power points in the sampling data set at 7.6 km seem 
to have imresolvable issues and have been omitted. As discussed in detail in the appendices, both 
spectroscopic methods derive emission indices from computations of mass fluxes of the chemical 
species of interest, a computation that required exhaust flow properties obtained from probe rake 
measurements. There is the potential for systematic errors in these conversions of number densities 
to emission indices, which might apply more or less equally over an entire altitude data set. It can 
be seen from an examination of the full IR diode laser absorption data set in Appendix D that in 
addition to the agreement with N’O measurement techniques, neither the H,0 nor the CO, emission 
indices arc systematically lower or higher than their theoretical values. The largest differences 
appear in the comparison of N0 2 values, to which we now- turn. 

Nitrogen species other than NO and NO, were not measured by the extractive gas: analysis 
(sampling) system. The probe sampling values for NO } concentrations, and thus F,I(N0 2 ), were 
derived from the difference of NO and NO, measurements with an implied assumption that no other 
nitrogen oxide .species contribute to the NO, measurement. On the other hand, measurements cf 
NOj concentrations by IR-TDL absorption are direct and no inferences arc required. However, as 
discussed in detail in Appendix C, the spectral region for this measurement has NO_, and H 2 () lines 
that are overlapped, such that a direct determination of the baseline (zero absorption) is not possible. 
Therefore, uncertainties may be large for IR measurements of N0 2 concentrations as well. IR 
measurements of NO, were made at only six engine test conditions, whereas the sampling measure- 
ments provide EI(NOj) values for almost the entire range of engine conditions . A quantity that is 
expected to more directly correlate with engine conditions is the ratio of EI(N0 2 ) to total EI(NO,). 
However, even for this ratio, variations between altitude data sets are such that at present the 
clearest trend is seen by plotting data for all altitudes and looking at the ratio as a function of T3. 
This plot, presented in Fig. 1 8, shows the N0 2 fractions for all but idle (lowest T3) conditions to be 
grouped around n canonical value of 0.1 , with NOj at idle being substantially iarger. Ratios between 
comparable cases for IR and sampling measurements range from 1 to 3. Ratios of absolute El(NO ; ) 
determined from IR and probe sampling measurements for the three lower altitudes with T3 = 733 
K are 1.34, 1.23, and 1 .03, agreeing well within the uncertainty limits. For two 15.2-km cases, the 
IR/sampling ratios of fsJ(NO ; ) arc 2 and 3, and for the 9.1-km idle case (where large differences 
exist for other species), the ratio is 2.5. 

An important aspect is the first-time demonstration of the viability of the IR-TDL technique 
to obtain high-precision, high-resolmion infrared absorption spectral measurements of an aircraft 
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engine exhaust over an extensive test period. Prior to the application of this technique to an earlier 
AEDC test, there was concern as to whether plume effects on the laser be»<n would allow 
measurements at all. Prior to the present test, it was still not certain that optical surfaces inside the 
test cell would remain clean and properly aligned. The harsh optical environment (including water 
spray, oil fumes, and high vibration and noise levels in the test cell) was not detrimental to 
providing quality measurements throughout the entire test period. 

The viability of the TDI, technique- for exhaust measurements demonstrated in this test 
program means that it can he further developed to give high-sensitivity measurements of 
concentrations of a wide variety of trace (and major) species, some of which have yet to be 
measured in exhausts by any technique. As implemented in this test, its ability to produce quan- 
titative values lor exhaust species concentrations and emission indices is dependent on 
characterization of the optical path plume properties provided by other independent measure- 
ments, such as the probe rake measurements of temperature, pressure, and velocity. Other imple- 
mentations of the TDL technique are possible, such as the use of multiple lines of sight (as wa* 
done with the IJV measurements of NO), multiple absorption lines of a single chemical species 
to determine temperature, and the possibility of choosing lines of trace and major species with 
the same dependency on temperature and pressure, so that trace species emission measurements 
could be self-referencing in the way that sampling measurements now determine El(NO) and 
EI(NO x ) by referencing to measured C0 7 concentrations. 

The OH laser absorption instrumentation performed well throughout the test period. The laser 
system was isolated from the direct noise and vibration of the test ceil by using a 30-m fiber optic 
cable to locate the laser in the basement of the test facility. The detectors and optics mounted inside 
the test cell survived the harsh test environment. The exhaust OH concentrations, however, were 
below the lowest detectable concentration for the exhaust conditions, 1 O' 2 cnr\ The laser was tuned 
to several OH lines throughout the test period, but spectral absorption features were never observed. 
Throughout the test period, the laser beam was intermittently redirected through a discharge lamp 
containing OH, which verified functionality of the system and ensured that the spectral scanning 
range was tuned to an OH line. Even though OH was not quantified, an upper limit concentration 
was determined. 

43 PARTICULATE AND AEROSOL CHARACTERIZATION 

The IJMR MASS was employed to characterize the aerosol emissions over the range of 
engine flight conditions described in this report. These data characterize the aerosol generated in 
the engine and that continued to exist at a distance of 12 cm downstream of the exit plane of the 
engine exhaust nozzle. This is the first database of its kind and provides particulate characterization 
(i.e., total concentration, aerosol, number- and mass-based emission indices, aerosol (soot) volume 
fraction, total aerosol and non-volatile aerosol size distributions and hydration (soluble mass 
fraction)) for such a cruise simulation. Using calculation techniques discussed in Appendix E, the 
mean aerosol number-based emission index calculated from all measurements is (2.2 ±0.7) x 10*'; 
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Che associated mean mass-based emission index is 0.012 ±0.001, (where the quoted uncertainty in 
both cases is the standard deviation of the mean). There is no strong correlation between emission 
index and thrust (measured as C0 2 concentration), altitude or combustor iniet temperature (T3) and 
pressure (P3). However, there appears to be a correlation between combustion efficiency and 
aerosol emission index observed from the ground idle data points, where elevated emission indices 
for aerosols, CO, and unburned hydrocarbons are recorded. The particle concentrations and 
emission indices measured in this test program were small compared with those measured on other 
aircraft engines using the UMR MASS. This correlates well with the extremely small smoke 
numbers reported in this study. Typical size distributions are log-normal and peak in the vicinity of 
20- to 40-nm diameter. The volatile component of the aerosol is measured to be small and in most 
cases can be considered negligible. Most of the ground testing to date using the MASS approach 
has shown similar log-normal distributions and small volatile components for a range of current 
engine configurations. A review of the soot inception literature indicates that aggregates of 
individual soot spheroids are usually found in probe-collected samples; it is generally assumed that 
the aggregates are formed during combustion as well. Electron micrographs of particles, which we 
have sampled by electroscopic precipitation from other combustion sources, do show aggregate 
character, with fractal dimension varying with combustion conditions. This characteristic of 
e,xhaust particle morphology warrants further study since it can influence the part iculatc’ s evolution 
in the atmosphere. A research project in which the MASS approach is combined with electron 
microscopy and optical analysis of soot particles extracted from the primary combustion zones of 
laboratory' flames would yield insight on the nature and influence of soot morphology in engine 
particulate emission characterization. 

A complete listing of the particulate characterization data on which this summary of results 
is based can be found in Appendix E. 

As part of the aerosol measurement component of this project, a C0 2 detector designed and con- 
structed by McDonnell Douglas was evaluated. The results were in excellent agreement with the 
AEDC gas sampling data, and the device will be used in conjunction with the UMR MASS for 
determination of emission indices directly in future ground tests. 

5.0 CONCLUSIONS 

This measurement program produced a comprehensive set of data that provides measure- 
ments of gaseous and particulate emissions at simulated sca-level-static and altitude conditions 
including cruise operation. This database can be used for atmospheric assessment studies, 
comparison of optical (1R and UV) non-intrusive diagnostics and conventional extractive gas 
sampling methods, and parametric measurements for SLS to altitude extrapolation sl^ liics. 

Data from the gas sampling system indicated trends with parametric variation in combustor 
inlet pressure and temperature as expected. Most variations observed in the spatial emissions 
distributions were in the outer regions of the exhaust flow and result from fan bypass air mixing 
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with the engine core air stream. The emission indices showed little spatial variation. The gas 
sampling data should serve as a good reference data set for evaluation of the innovative optical 
techniques, also applied as a part of this test program. 

The optical linoof-sighc non-intrusive techniques provided chemical species concentration 
measurements independent of the probe sampling measurements. Overall agreement of the data 
collected from these independent measurement systems throughout the test period substantiates 
accuracy within respective uncertainty limits. Additionally, highly repeatable measurements during 
steady-state test conditions from all measurement systems provide evidence of engine operation 
stability, spatial uniformity, and radial symmetry of the exhaust. 

The IR-TDL system performed reliably throughout the test period, providing a large body of 
data on NO, NCX CCX and H : ,0 concentrations. The method used £o interpret observed infrared 
absorption optical depths (using a model of the single nonuniform line of sight) produced trace 
species emission index values that agree with the other measurement techniques to within expected 
uncertainty limits, as well as major species emission indices that agree with theoretical values with- 
in similar limits. 

The NO UV resonance absorption technique with multiple LOS measurements across the 
exhaust flow field provided radial profiles of NO concentrations slightly lower but in good agree- 
ment with radial profiles measured using single-point probe sampling. Agreement was achieved 
even in the surrounding gas-recirculation region. Emission indices agreed over all test conditions, 
well within uncertainty limits for respective measurement techniques. 

Attempts at direct measurements of OH concentrations were unsuccessful due to nozzle exit 
OH concentrations below 10 12 cnr\ the detectable limit for the laser system as configured for these 
test programs. However, this measurement establishes an upper bound on OH concentration for 
these, tests which is important information for ncar-ficld plume interaction studies. Also, since 
accurate OH concentration prediction models are not available, these results provide design criteria 
for developing accurate chemical mechanisms for OH production in turbine engine exhausts and 
designing experimental simulations of these environments. 

It is important to note that each optical measurement syst. m produced precise and repeatable 
spectral data; individually predicted overall consistent trends; and agreed wifi, in stated uncertain- 
ties of each instrument and probe sampling technique. Although the uncertainties were reasonably 
low, the greatest uncertainty contribution was attributed to the characterization of the optical path 
properties required to deduce densities and emission indices for respective species. 

The following major conclusions can lx deduced from the aerosol characterization data. The 
engine studied in this* test is representative of most of the engines characterized in ground tests, to 
date, with this approach. Namely, the nature of the size distributions, being log-normal, is typical, 
ns is the mean size of the distribution and its width. The third parameter characterizing a log-normal 
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distribution, the total number of particles, is the lowest measured for any engine characterized to 
date. The similar result for smoke number measured for this engine suggests a correlation exists 
between smoke number and particle emission index, at least in the case where these numbers arc 
small. In addition, the volatile component of the particles is small. This is also consistent for other 
ground-based measurements where the sampling probe was located less than 1 m from the exit 
plane of the exhaust nozzle. 

The particle emission indices were not determined to be strongly dependent on either fuel/air 
ratio (as measured by CO ; concentration), altitude, or combustor inlet pressure or temperature. 
Based on the limited data set: for the ground idle condition, where the engine combustion efficiency 
is at its lowest, there appears to be an inverse relationship between the particle emission index and 
combustion efficiency. The implications of this result with respect to our understanding of particle 
production as a function of engine performance warrant further investigation. The apparent absence 
of aggregates in the particulate emissions detected with the MASS raises interesting and critical 
questions with respcci to our understanding of the soot formation process and/or the ability of the 
MASS methodology to account for their presence. 

Overall, this emission test program exceeded expectations in the quantity of data and engine 
test conditions achieved for the limited test period. The test demonstrated the viability of several 
diagnostic techniques to characterize critical emissions parameters of interest and thus provid . 
comprehensive set of turbine emissions data to the scientific community for aviation impact 
assessments on atmospheric chemistry. Further measurement programs will be required to develop 
a comprehensive database spanning all classes of engines in current and future aviation fleets. 
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Figure 3. Horizontal COj profile for T3 = 733 K at 9.1 -km altitude. 
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Figure 5. Horizontal EI(NO) profile for T3 - 733 K af 9.1-km altitude. 
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Figure 12. Horizontal bar EI(THC) summary for all test conditions. 


29 













- Sampling 


AEDC-7R-9&-3 


c 

o 

•£ 

O 

</> 

< 

4) 

2 

m 

c 

8 

£ 

?D 


i 0 -i 


o 

to 


V It 

I ^ 

\ 

hWWM 

V. 

• ^ 

' V 
> ^ 

(■4-ic-V-H 




L 


o 


<M 


-L 

o 


U R 


( Z 0H sb) xepuj U0!5*IU13 ON 





AEOC-Tfi-96-3 



34 


500 600 700 800 

Combustor Inlet Temperature (T3), K 

Figure J7. Comparisons of EHNO) at 15.2-km altitude. 
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Table 1. AEDC Fuel Analysis Report 


Laboratory 

Laboratory Test 

— — 

Laboratory Results 

Analysis 

Method 

Pretest 

Midtest 

Posttcst 

ASTM (1) D2382 

Heat of Combustion, btu/lbm Gross 

19814.6 

19818.0 

19820.8 

ASTM D2382 

Heat of Combustion, btu/lbm Net Min. 

18544.7 

18547.2 

18550.0 

ASTM D3701 

Hydrogen Content, percent by Wt. 

13.92 

13.93 

13.93 

ASTM D3828 

Flash Point, °F 

144 

143 

145 

ASTM D4052 

Specific Gravity ( 1 5C/4C), Min. -Max 
@ 60°F/4°C 

08147 

0.8147 

0.8148 


@100 o F/4°C 

0.7985 

0.7986 

0.7987 

ASTM D445 

Viscosity @0°F, cs 

mm 


5.730 


70"F. cs 

B|| 


2.165 


150°F, cs 

■ESI 

1.115 

1.117 

TO< 2) 42B-l-l 

Free Water, PPM, max 

0 

0 

0 


American Society for Testing and Materials 
Air Force Technical Order 
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Table 2. USAF Aerospace Fuels Laboratory Fuel Analysis Report 


Laboratory Analysis Method 


ASTM* D4176 


ASTM D3242 


AS'FM DI319 


ASTM D1319 


ASTM D3227 


ASTM D4294 



ASTM D3338 


ASTM D3343 


ASTM D1322 


ASTM D 130 


Laboratory Test 


Workmanship 


Total Acid Number, mg KOH/g 


Aromatics, percent vol 


Olefins, percent vol 


Mercaptan Sulfur, percent wt 


Total Sulfur, percent wt 


Distillation IBP, °C 
10 percent 
20 percent 
50 percent 
90 percent 
EP 

Residue, percent vol 
Loss, percent vol 


Flash Point, e C 


Gravity, API 


Freezing Point, °C 


Viscosity <s ; -20 9 C, cs 


Heat of Combustion, MJikg 


Hydrogen Content, percent wt 


Smoke Point, mm 


Copper Strip Corrosion 


Laboratory Results 
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Table 2. Concluded 


Laboratory Analysis Method 

Laboratory Test 

Laboratory Results 

ASTM D3241 

Thermal Stability 
Tube Deposit visual 
Change in Press., nun of Hg 

i 

0 

ASTMD381 

Existent Gum, mg/l 00ml 

0 

ASTM D5452 

Paniculate Matter, mg/l 

0.1 

SPEC 

Filtration Time, min 

10 

AS'Cvl D1094 

Water Reaction 

2 

ASTM D5006 

FSII, percent vol 

0.16 

ASTM D2624 

Conductivity, pS/m 

244 

ASTM D5001 

Rall-On-Cylinder, mm 

0.56 


Trace Metal Results 
! Cu 

| Fe 

Zn 

m 


“Arnericiui Sin;it:l> for Testing and Materials 


.18 
































Table 3. Engine and Flight Conditions 
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Table 4. Measurement Systems Used to Characterize Exhaust Emissions 


Measurement 

Measurement 

Principal 

System 

Quantity 

Agency 

Optica! Non-Intrusive 



IRTunable-Diode Laser Absorption 

NO, NO*. C0 2 , H 2 0 

Aerodyne Research. Inc. 

NO-UV Resonance Absorption 

NO 

Sverdrup/AEDC 

NO-UV Continuum Absorption 

NO 


UV-Lascr Absorption 

OH 





Extractive Sampling 



ARP1 256* 

C0,C0 2 . 0 2 . NO* NO x .THC 

Sverdrup/ AH DC" 

Smoke Number Meter* 

Plume Visibility (Soot) 


Hygrometer 

Inlet Air Humidity 


j 

Mobile Aerosol Sampling I 

Particulate/ Aerosol: Total Mass, 

University of Missouri- 

System (MASS) 

S J2-C Distribution, Hydration Properties, 
Emission Indices 

Rolla 

Extractive IR-Absoiption 

C0 2 



McDonnell Douglas 




Probes 

Sample Probe 

Exhaust Sample Extraction 

Svcrrirup/AEDC 

Local Mach Number. Flow Angularity, 


Cone Static, Pilot Pressure, and 
Total Temperature Probes 

Static Temperature, nnd Static Pressure 

S verdi up/AEDC 





4 Aerospatr Recommended Practice (AKIM Y?St B ’'PiuccJimr foi the Continuous Sampling and Mcastiremcnt of Gaseous 
Emissions fronitia* Turbine Entities'* 


41 













Table 5. Summary of Exhaust Emissions Characterization 
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APPENDIX A 
PROBE RAKE SYSTEM 

R.S. Hiers, Jr. 

A cruciform configuration rake was located just downstream of the engine nozzle exit plane 
in the test cell diffuser/hcat shield to provide extractive emission and jet plume aerodynamic da':a. 
A photograph illustrating the rake installed behind the engine is shown in Fig. A-l. 

A- 1.0 SAMPLE EXTRACTION PROBES 

The primary function of the rake was to provide a well-resolved, spatial distribution of 
extractive gas samples to a suite of SAE qualified/standard emission instrumentation for 
characterizing engine emissions. For this purpose, the extractive gas sample probes were located on 
nominal 5.08-cm centers on both the vertical and horizontal members of the cruciform rake, 12.7 
c-m downstream of the nozzle exit plane. For completeness, extractive probes were positioned in the 
hot core flow, the bypass flow, and also in the test cell cooling flow external to the engine flow path. 

A sketch of a gas sample extraction probe is shown in Fig. A-2, illustrating the probe's sharp- 
edged inlet orifice and nominal 10 to I area expansion just downstream of the inlet. These features 
eliminate mutual mass flow interference between interconnected (ganged) probes when suitable 
sample system back pressures are maintained on the extractive probes. 

The sample transport lines from each extractive probe were routed individually to a valving 
system which permitted both individual gas samples or collective (ganged) gas samples from 
desired groupings of individual extractive probes to be obtained. 

The sample transport lines were heated with steam tracing and insulated to maintain a 
temperature of approximately 422 K, from the location where the sample transport line emerged 
from the engine exhaust heated portion of the rake structure to the valving system, to prevent 
condensation. Of course, on the rake structure influenced by the hot engine flow, extractive sample 
lines arc typically well above 422 K. 

The valvir.g system and the sample transfer lines connected directly to the emission analyzers 
were also maintained about 422 K with electrical heat tracing and insulation. 
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A-2.0 AERODYNAMIC PROBES 

Mach/flow angularity and stagnation temperature probes were also employed on the rake for 
providing local pointwise static pressure, static temperature, density, and velocity along the line-of- 
sight (LOS) of both the UV and the !R optical instrumentation. The optical instrumentation provid- 
er! "non -intrusive" measurements for comparison and validation with extractive probe data (NO 
concentration primarily). 

In the past, computational model results were used to provide the local thermodynamic 
conditions in the LOS of optica! instrumentation. Local gas properties are used for molecular 
pressure broadening corrections and energy level populations, and for computing relative density 
(e.g., ppm) from the absolute densities directly. 

In the present work, the thermodynamic parameters were measured to eliminate the 
uncertainty inherent in using model results and consequently to provide better assessment of the 
fundamental accuracy of these optica! techniques. A secondary by-product of these data is the 
ability to validatc/assess the engine plume models by comparison with rake data. 

Seven conventional Mach/flow angularity (MFA) probes and 10 stagnation temperature 
probes were mounted on the horizontal member of the rake, interspersed between the gas extraction 
probes. Design details of the MFA and stagnation temperature probe are shown in Figs. A-3 and A- 
4. Tire MFA probes and the total temperature probes were all located on 10. 16-cm centers. The 
physical, relationship and spacing (2.54-ctn probe-to-probe) between the extraction, MFA and total 
temperature probes is illustrated in Fig. A-5. 

A -3.0 DATA REDUCTION 

The basic data reduction procedure for inferring the free-stream velocity and static 
thermodynamic flow properties from the combination of MFA, stagnation temperature, and 
extraction measurements is illustrated schematically in Fig. A-6. 

Briefly a computational fluid dynamic (CFD) relationship was developed relating the ratio of 
average cone surface pressure to cone pitot pressure of the MFA probe to the frcc-strcam Mach 
number. The relationship is strongly dependent on MFA probe configuration and as evident from 
Fig. A-5, Mb' A probe geometry was carefully controlled during manufacture and testing. The 
relationship of the cone surface to cone pitot pressure ratio to Mach number is much less sensitive 
to composition. Nevertheless, the composition determined from the extractive probes was used to 
precisely describe the gas internal energy characteristics for the correlation of MFA probe pressure 
data with Mach number. 
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Once the local Mach number just upstream of an MFA probe location was inferred from 
MFA probes as described, the local static pressure was inferred from the gas dynamic relationship 
between the measured cone pitot pressure and the now known Mach number. 

Static gas temperature upstream of MFA probes was similarly inferred from the gas dynamic 
relationship between the now known Mach number and the measured total temperature (again 
accounting for the minor effect of gas composition). 

Unfortunately, since the Mach number and stagnation temperature were not measured 
coincidentally in space, the total temperature that would actually exist at an MFA probe required 
estimation by interpolating from measured total temperatures on either side of the desired probe 
location. 

Once temperature and pressure (along with composition) were determined, all the remaining 
thermodynamic properties, e.g., density, were unambiguously computed. Local stream velocity was 
also determined from Mach number and static temperature, completing the definition of the local 
state of the flow field at a discrete point in the LOS of the optical instrumentation. 

A-4.0 RAKE DATA UNCERTAINTY 

Typical uncertainties in the measured probe temperatures and pressures are of the order of 
0.75 percent and 0.25 percent, respectively. The maximum uncertainty in any of the derived 
thermodynamic parameters is about 2 percent as a result of propagation of measurement 
uncertainties through the data reduct ion algorithms. The uncertainties arc primarily bias errors. The 
use of a CFO-generated “calibration" for the Mach number probes, as previously described, 
introduces an additional source of uncertainty compared with the use of an actual aerodynamic 
calibration for the MFA probes. However, this technique has been previously used and validated at 
AEDC for similar conditions where "CFD calibrations" were compared with actual aerodynamic 
MFA probe calibrations. For larger MFA probes, typified by the present design where the externa! 
aerodynamic geometry can be precisely controlled during manufacture and test, CFD calibrations 
have been found to be essentially equivalent to aerodynamic calibrations for the spherically blunted 
class of conical MFA probes. Also in previous work, it was found that for other MFA probe 
geometries or much smaller probes, aerodynamic calibrations were required because smaller 
probes’ geometries cannot be rcpeatably manufactured. 

Since the large-scale class of MFA probes was used for these tests, the uncertainty introduced 
by the use of CFD calibrations is believed small and is conservatively estimated to result in a total 
uncertainty in any of the rake aerodynamic properties of no more than 5 percent, based on the 
historical experience described above. 
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A validation of (he quality of the rake data is illustrated by comparing the mass flow 
integrated from rakc-derived velocity and densities with the independent measurements of total air 
and fuel flow provided by the test facility instrumentation. 

In spite of the interpolation required of the relatively sparse probe data and the extrapolation 
required in the region of high gradient near the edge of the engine nozzle, the integrated mass flow 
is clustered within ±15 percent of the total mass flow to the engine over the range of power and 
altitudes tested. 

Because of the large extraneous uncertainties introduced into these comparisons by the 
interpolations and extrapolations required for the integration, the above comparisons suggest that 
the estimated uncertainty in individual rake properties of less than 5 percent is realistic. 

Of course, if this becomes an issue in analysis of LOS optical data, the MFA probe may yet 
be calibrated aerodynamically and the uncertainty more accurately determined. The best use of an 
aerodynamic calibration would be to re-reduce the aerodynamic data with a concurrent lowering of 
the residual uncertainty. This remains an open, attractive option. 

A-5.0 RAKE AERODYNAMIC DATA 

Rake aerodynamic data for all test conditions are summarized in Table A- 1 . The location of 
the probes in the engine exhaust plume is relative to the nozzle exhaust coordinate system, shown 
in Fig. A-7. 
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Figure A-5. Typical probe installation on horizontal rake. 
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Figure A-7. Gas sampling probe locations near the nozzle exit. 
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Table A-l. Continued 
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APPENDIX B 

EXTRACTIVE GAS ANALYSIS SYSTEM 


V. A. Znccardi, D. G. Gardner. R. P. Howard, D. J. Gonnion, 

P. A. Jalbert, J. A. Ready, R. S. Hicrs, Jr. 

B-1.0 INTRODUCTION 

Conventional extractive gas sampling techniques arc used for engine emissions 
certification and thus are accepted as standard methodology (Refs. B-J, B-2, B-3) for character- 
ization or gaseous constituents. The system used for this NASA test program is similar to a sys- 
tem that would be used for commercial aircraft engine emissions certification testing and 
therefore should provide reliable data for this study and comparisons to optical techniques also 
used for species concentration measurements. 

B-2.0 INSTRUMENTATION 

Concentrations of gaseous species contained in the engine exhaust were measured using a 
suite of gas analyzers mounted outside the test cell. Constituents measured were nitric oxide (NO), 
oxides of nitrogen (NO,), carbon monoxide (CO), carbon dioxide (C0 2 ). oxygen (0,), and total 
unbumed hydrocarbons (THC). Nitrogen dioxide (N0 2 ) concentrations reported in the main body 
were inferred from the difference of NO, and NO concentration measurements, assuming other 
oxides of nitrogen were negligible. In addition to gaseous components, measurements of black 
smoke were measured and reported as smoke number. 

A cruciform multi-point rake, described in Appendix A, was used to obtain exhaust 
samples at individual and combined probe groupings that provided radial profiles of species 
concentrations as well as spatial averages. An illustration of the rake mounted just inside the test 
cell diffuser is presented in Fig. 2 of the main text. The cruciform rake was mounted inside the 
diffuser with the tips of the gas sampling probes located about 12 cm downstream of the nozzle 
exit plane. The rake system and design of the probes used for sampling are described in Appendix 
A. Probe locations with respect to the diffuser and nozzle exit area (approximating exhaust flow- 
field area) arc presented in Fig. B- 1 , a view looking downstream. Probe numbers in Fig. B-l were 
used to designate probes sampled during data acquisition. The rake was made oi two bars, one 
mounted vertically and containing only gas sampling probes, and the other mounted horizontally 
with sampling probes interspersed among Mach/flow angularity (MFA) and total temperature 
probes. The center of the rake, and thus the horizontal bar, was 2.54 cm above the centerline of 
the nozzle. The sampling probes were equally spaced at 5.08-cm centers, except lor a 10.16-cm 
distance between the center-most probes (numbers 10 and 28) on the vertical bar. There were six 
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probes on each radial arm of the horizontal bar and lower radial arm of the vertical bar contained 
within the exhaust flow field. Only five probes of the upper radial arm of the vertical bar were 
within the exhaust flow field. 

Probes 7-9 and 25-27 were used for measurements of species concentrations in the region 
surrounding the exhaust flow field where exhaust recirculation gases could interfere with optical 
technique measurements. Probes 7 and 25 were part of the horizontal bar assembly. Samples at the 
other locations (outside the exhaust flow field) were extracted through stainless-steel tubing 
mounted either to the rake braces or other hardware in the test cell. 

Gas samples collected at the probes were transported through stainless-steel lines to a gas 
distribution system. The gas sampling lines (one for each probe) were heated, insulated, and 
wrapped to maintain the temperature at 422 K. Inside the cel!, the sample lines were heated with 
steam lines. Outside the cell, the lines were electrically heated using heat tape. The total length of 
the sampling lines from the rake probes to the gas distribution system was about 1 5 m. 

The gas distribution system had four manifolds, one for each radial arm of the horizontal and 
vertical bars. These manifolds were designated as 90°, 1 80°, 270°. and 360° as illustrated in Fig. B- 
2. Each sample line contained a valve that either directed the sample into a manifold, or into a 
sample bypass line where the sample was discharged and discarded. The output of each manifold 
was directed to either the AEDC gas sampling instrumentation (discussed here), to the University 
of Missouri-Rolla Mobile Aerosol Sampling System (1JMR-MASS) discussed in Appendix E, or to 
the sample discharge line. Simultaneous measurements with independent samples were achieved 
for gaseous and particulate characterization with the restriction that AEDC and UMR would not si- 
multaneously select the same radial arm (manifold). As stated earlier, this arrange ment allowed sin- 
gle-point and multiple-point averages for each measurement system simultaneously. The entire gas 
distribution system was heated and maintained at 422 K. 

Gas samples were transported a short distance from the gas distribution system to the suite 
of gas analyzers through a single heated line. The gas sample was branched into several analyzers 
simultaneously for test conditions that provided enough flow. Otherwise, measurements were made 
sequentially by flowing the gas sample into subgroups or individual analyzers. 

Manufacturers, instrument model numbers, and span gas concentrations are given in Table 
B-l for the analyzers used for measurements of gaseous species in the exhaust samples during this 
test program. Non-dispcrsive infrared (NDIR) analyzers (dry) were used for measurements of CO 
and C0 2 concentrations. An analyzer (dry) using a paramagnetic detector (PMD) was used for 
measurements of 0 2 concentrations. Chemiluminescence (CLD) analyzers (dry) were used for mea- 
surements of NO and NO*, with the sample passed through a catalytic converter for NO* measure- 
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merits. One of the- two chemiluminescence analyzers failed early in the test period requiring that 
both NO and NO, be measured sequentially on the one working analyzer for each data point. Total 
unbumed hydrocarbons (THC) were measured (wet) using a heated flame ionization detector (FID). 
Except for THC, all measurements were performed on dry samples. Dry indicates that the samples 
were passed through a gas sample dryer (Perron Pure Dryer-Model PD- 1000) before entering the 
analyzer. A cold mirror dew-point hygrometer was used for real-time monitoring of dryer 
efficiency and provided humidity correction factors for carbon balance calculations. 

National Institute of Standards and Technology (NTST)-traceable calibration gases were used 
for routine verification and adjustment of the analyzer response for each instrument. Gaseous 
nitrogen (GN 2 ) was used to establish the zero or non-response of the analyzer. The primary 
measurement ranges and calibration gas concentrations used to adjust the analyzer response are 
presented in Tabic B- 1 . 


The black smoke meter was a standard Aircraft Gas Turbine Engine Exhaust Smoke 
Measurement system meeting the requirements of the SAE Aerospace Recommended Practice 
(ARP1 1 79B, Ref. B-2). The schematic in Fig. B-2 shows the system in relation to the. gas analyzers. 
Gas samples collected by the probe rake were metered through a filter paper. The change in reflec- 
tance of the paper due to the smoke particles was related to a smoke number (SN). 

B-3.0 DATA ACQUISITION 


Maximization of the number of engine test conditions without sacrificing data quality was 
achieved through diligent real-time examination of the gas sampling data. Strategic data points 
planned early in the test period provided assurance of engine and gas sampling system stability arid 
radial symmetry in the exhaust flow field. Radial symmetry was observed independently on the 
vertical and horizontal sampling bars during detailed mapping of the engine core, fan bypass, and 
intermediate mixed regions at steady-state engine test condition. This level of detail was repeated 
at priority engine test conditions at each of the several altitudes. Otherwise, the horizontal bar was 
devoted to gas sampling data, while the vertical bar was devoted to particulate characterization 
measurements, described in Appendix E. 

The ambitious test matrix planned for a single lest period required automation of the 
electronic data acquisition system to provide real-time concentrations for step-by-step decisions 
throughout the test program. Automation was achieved using LabVIEW* to control data acquisition 
of each analyzer system. Raw data from the gas analyzers were saved to the hard drive of the data 
acquisition personal computer (PC) and to the hard drive of the data reduction PC. via an RS232 
line. The separate computer allowed data reduction without interruption of LnbVIEW during 
analyzer monitoring and data acquisition. The Lab VIEW data file was accessed directly by the PC 
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supporting Microsoft ExceP for immediate data reduction (Fig. B-3). An Excel spreadsheet was 
used for calculation of the emissions parameters. In addition to reducing raw data voltages to 
respective species concentrations, preliminary carbon balance calculations were performed for each 
gas sampling data point. The Excel-calculated emissions parameters were printed immediately and 
a copy provided for real time data assessment as explained earlier. The raw voltages and standard 
deviations of the measured species data set were included in the report for assessment of the 
analyzer calibrations and sample stability, respectively. 

B-4.0 RESULTS 

All emissions data arc presented in Table B-2 which also contains the performance data 
required for analysis and designation of engine conditions. The data are presented in order of 
acquisition during the test period. The first two columns show the sample probe configuration 
related to the gas sampling data point (DP) number in column three. The radial arms of the cruci- 
form rake are referred to as the left (270°) and right (90*) horizontal arms, and the upper (360°) and 
lower (1 80°) vertical arms; abbreviated LH, RH, UV, and LV, respectively. An “x” in the table des- 
ignates active sampling from the respective probe position (looking downstream), and “o'* indicates 
the probe was not sampled for that data point. The brackets [ ) enclose probe positions within the 
nozzle exit region. Thus, the first data point in the table indicates that all left horizontal (LH) sam- 
pling probes within the exit area were sampled simultaneously for DP6. Probe designators outside 
the bracket in the first column correspond to probe positions outside the engine exit, but in the same 
horizontal plane. The physical positions of these probes relative to the nozzle exit plane and 
entrance lip of the diffuser are given in Fig. B- 1 . Vertical probes outside the engine exit diameter 
were not sampled and thus probe designators were not necessary outside the vertical arm brackets 
in column two. 

The clock time for each data point is given in column four. The altitude (ALT), Mach num- 
ber (MACH), combustor inlet total temperature (T3), and combustor inlet total pressure (P3) re- 
late the data point to the steady-state engine condition. Species concentrations are reported as 
fractional densities by volume with major species in percent and minor species in parts per mil- 
lion volume (ppmv). The species concentrations [OJm ( [C0 2 ]. [CO], [NO,], and [NO] 

were measured. The quantities [0 ? )c. fH 3 OJ, and f N^] were calculated from the method of carbon 
balance described in Ref. B-3. Since oxygen values were both measured and calculated, these 
values were denoted using "m" and “c'\ respectively. The measured (QJm data agree well with 
the carbon balance calculations. Although agreement docs not guarantee accurate data, signifi- 
cant deviation would have indicated a problem for the respective data point, The emission index 
(El) for each measured minor species was calculated according to Ref, B-3 with EI(NO) reported 
as t v *c N0 3 equivalent, as defined in Ref B-3. The fuel-to-air (fud/air) ratio and combustion ef- 
ficiency (percent) were determined by the method of carbon balance for each gas sampling data 
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point. Fuel/air ratios were determined independently from engine performance data for both the 
engine hot core and mixed core and fan bypass air and are presented in the table for comparisons. 
Finally, smoke number, measured at only three engine test conditions, are given in the last col- 
umn. All smoke numbers were less than 1 .5, which is within the zero scatter of this technique; 
thus indicating that smoke was not measurable. 

REFERENCES 

B-f. Procedure for the Continuous Sampling and Measurement of Gaseous Emissions from 
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Figure B-l. Gas sampling probe locations near 
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Table B-l. Summary of Emissions System Analyzers and Span Gases 


TYPE 

MODEL 

RANGE 

CALIBRATION 
SPAN HAS 

CO-NDIR 

California Analytical, Inc. 
Model ZRH 

0- 1000 ppm v 

826pprnv 

COj -NDIR 

California Analytical, Inc. 
Model ZRH 

0-5 percent 

3.99 percent 

NO-CLD 

Thermal Electron Corp. 
Model 10 

0-25 ppmv 
0-250 ppmv 

203 ppmv 
202 ppmv 

NO x -CLD 

Thermal Electron Corp. 
Model 10 

0-25 ppmv 
0-250 ppmv 

203 ppmv 
202 ppmv 

THC-FTD 

Beckman 
Mode) 402 

0-100 ppm vC 
0- 1,000 ppmvC 

89.7 ppmvC 
885 ppmvC 

0 3 - Paramagnetic j 

California Analytical, Inc. 
Model 5000 

25 percent 

20.8 percent 
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Tabic B-2. Gaseous Emissions and Smoke Summary 
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APPENDIX C 

NO-UV RESONANCE AND CONTINUUM ABSORPTION SYSTEMS 
R. P. Howard and A. G. Davis 
01.0 INTRODUCTION 

Spectral interrogation of nitric oxide (NO) resonance electronic (A«-X) transitions using 
ultraviolet (UV) absorption allows determination of NO density and under appropriate 
conditions, thermodynamic state. NO-UV absorption techniques, developed in the 1970’s (Ref. 
C-l) for measurements of low NO concentrations in gas turbine combustors and engines, were 
modified and enhanced in response to new requirements for non-intrusive measurements in high 
enthalpy test facility flows (Refs. C-2 and 3), exhausts of liquid-propellant rockets 
(unpublished), and for automobile exhausts (Ref. C-4). The hardware was made robust and more 
reliable using spectrometer focal plane array detectors, fiber optics, and pulse-start power 
supplies for the resonance lamps. The line-by-line radiative transfer model, previously limited to 
the NO gamma (0,0) band, was extended to include the (GJ) band and thus allowed simultaneous 
determination of NO number density and ground electronic state-vibrational population 
distribution (reported as ground state- vibrational temperature) in many of the aforementioned 
applications. Renewed interest in effects of turbine-powered aircraft on atmospheric chemistry 
has again prompted exhaust characterization studies and new requirements for measurements of 
NO in turbine engine exhaust. 

For this test program, NO-UV resonance and continuum absorption techniques offered in- 
situ non-intrusive optical measurements of NO number density in the engine exhaust over the 
full set of operation conditions. The NO-UV line-by-line radiative transfer model coupled to a 
radial inversion scheme was used to determine radial NO number density profiles from spectral 
transmission measurements through several chords of the turbine exhaust. 

These measurements were sponsored by NASA under the Atmospheric Effects of Aviation 
Project (AEAP) as a part of a parametric study of exhaust emissions of a gas turbine engine 
operated over a wide range of engine conditions at each of several simulated altitudes. NO-UV 
resonance absorption and continuum absorption measurements were made at each steady-state 
condition. Described here arc the measurements, data reduction to transmittances, application of 
the NO radiative transfer model and inversion scheme for determination of NO number density 
profiles, and determination of an NO emission index for each steady-state test condition. 
Descriptions of the engine and facility parameters for each steady-state flight condition are 
detailed in Section 2.0 of the main text of this report. The resonance and continuum absorption 
system components were similar except for radiation sources- The system and technique 
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similarities make it convenient to describe system components, data acquisition, and data 
analysis for both systems simultaneously while pointing out differences only as necessary. 

C-2.0 INSTRUMENTATION 

The resonance and continuum system instrument components differed only in lamp types. 
Primary components per system included a lamp, a lamp shutter with an elec tronic controller, 
two fiber optics cables, collimating and focusing optics, a grating spectrometer, an optical 
multichannel analyzer (OMA*) detector system with a computerized controller, and a personal 
computer for data acquisition and storage. The systems shared stepper motor driven slide tables 
fortraversing the optical beams across the exhaust flow field. 

The resonance lamp, a source of NO resonant radiation, was a DC-excited capillary 
discharge tube operated at 12 ma using a 5,000- V power supply with a 220-kohm ballast 
resistance. A gas mixture of argon:nitrogen:oxyg.en, at a 12:3:1 ratio by volume, was flowed 
through the lamp tube while maintaining a static pressure of 10 torr. The capillary discharge tube 
was water cooled to tap water temperatures. Resonance lamp intensity characterization is 
discussed in Section 3.0 below. The continuum system lamp was a constant voltage 100-W, 
water-cooled deuterium lamp. 

Each detector system consisted of a 0.32-m Czerny Turner spectrometer with a variable- 
width entrance slit, a 2,400 gnxjve/mm grating fused in first order), and an EG&G Princeton 
Applied Research Model 1421.-B silicon-diode array detector. The detector was UV-enhanced 
and intensified over the full 1,024 pixel elements. The detector and detector controller (model 
1461) arc referred to as the OMA. The spectral range was 214 to 245 nm with a 0.029-nm per 
detector-element spectral dispersion and 0.16-nm spectral resolution (full-width-half-maximum). 
The instrument response function for a single array element, g(\) in Eq. (G-l) of Section C-3.0, 
is presented in Fig. C-l for the resonance absorption system. 

The spectral region reported for the continuum absorption measurement deserves 
explanation. First, this measurement system was developed and installed with intentions of 
measuring OH resonance absorption near 300 nm. After successfully incorporating a UV 
spectral- scanning laser for OH absorption and having failed to detect OH absorption during an 
engine checkout run, this system was reconfigured to allow NO continuum absorption 
measurements. The combination of system spectral throughput (primarily the long fiber cables) 
and the deuterium lamp characteristics resulted in detector saturation in the region of the (0,1) 
band for detector integration times required for sufficient signal-to-noise ratio at the (0,0) band. 
Therefore, the less important (0,1) spectral region was moved off the detector. Also, there was 
insufficient lamp signal for meaningful measurements at wavelengths below the (0,0) band. The 
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useable spectral range was about 222 to 228 nm with a 0.029-nm per detector-clement spectral 
dispersion and 0.21-nm spectral resolution (full-width-half-maxi mum). The instrument response 
function for one array element of the continuum system, corresponding to g (X) in Eq. (C-l), is 
also presented in Fig. C-l . 

C-3.0 RADIAL INVERSION-RADIATIVE TRANSFER MODEL AND 
RELATED PARAMETERS 

A computer model based on the theoretical line-by-line radiative transfer for the NO 
gamma (0,0) band was used to determine NO number density from spectral absorption 
measurements. The theoretical physics of NO molecular absorption, details of the model, and 
extension of the model from earlier work can be found in the references (Refs. C-l, 2, 5> 6). A 
brief overview is presented here to describe and document model parameters used for analysis of 
the present data. Basically, the fractional lamp radiation within the wavelength interval, AX, 
transmitted through a medium of length, L, is given by the equation of radiative transfer, 

w \ -Y 

M Zj<'c + '»;<*)>*«'»■ <c ’ l) 

j 

where 

g(X) = spectral instrument response function, 

/r^ f (X) = molecular absorption transition (line), 

/>/. (X) = relative intensity of a resonance lamp transition, and 

/ t = Lamp radiation factor, modeled as continuum radiation. 

The interval AX is limited by the spectral extent of g(X), the spectral response for an individual 
element of the linear array detector. Absorption in the test media is modeled by molecular 
absorption transitions, A>,(X), which arc described by Voigt line-profiles that allow for pressure or 
collision broadening. The summation over i includes all transitions of the NO gamma (0,0) or 
(0,1) band wholly or partially falling within AX. The absorption coefficient is a function of the 
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static temperature, static pressure, directed velocity line center shifts (Doppler shifts), line 
center pressure shift parameter, broadening parameter (discussed tater), and N'O number 
density along the measurement path, L. For a single homogeneous path, the flow-field 
parameters are input to the model, and the NO number density is varied until the calculated 
(Eq. (C-l)) and corresponding measured transmittances agree. Calculations were compared to 
peak absoiption at the second bandhead, 226.23 nm and 236.3 nni, respectively, for the (0,0) 
and (0,1) bands using the resonance lamp and 226.18 nm for the (0,0) band of the continuum 
absorption technique. The (0,1) band was not measured with the continuum absorption system, 
as discussed in Section C-2.0. 

Absorption media line center shifts with respect to the fixed lamp radiation line positions 
are important to the NO densities derived from the model. Two shift mechanisms are modeled in 
the code. Shifts from the Doppler effect were negligible for these measurements due to small 
exhaust flow velocity components in the direction of the source beam. Collision-induced shifts 
were modeled according to the collisional shift parameter, 5 = -0,18 V P(295/T)" 5 \ as quantified 
in Ref. C-l. Here P is static pressure in atmospheres and T is static temperature, in Kelvin. 

The resonance lamp radiates several lines within the spectral interval, AX, of Eq. (CM). 
Line intensities, !- A ] (X), are modeled by Doppler profiles broadened (o 950 K. The relative 
lamp line intensities were measured and characterized according to the upper-state population 
distribution. Fig. C-2. Results for the (0,0) and (0.1) bands were similar, as expected, since the 
upper states are common to both bands. The lamp emits low-level radiation (l t ) in addition to 
gamma band radiation. This radiation was not directly quantifiable from lamp characterization 
measurements, but was modeled as constant continuum radiation, I k ., over the spectral range, AX, 
with excellent results. 

The lamp radiation for the continuum absorption technique was modeled as a simple 
spectral constant, achieved in Eq. (C-l) by setting l- K * (X) to zero and 1 to some arbitrary 
constant value. 

For the resonance technique, the molecular broadening parameter, a' = C*l’/T n , and I. 
have been treated as instrument-dependent code calibration factors (Ref. C-2.) selected to 
minimize differences in calculated and measured transmittances at controlled laboratory 
conditions. Here, C is referred to as the broadening constant and n the temperature-dependent 
exponent. The broadening parameter, a' = 31800*P/T'\ was reported in 198) (Ref. C-5) but did 
not include a low-level lump radiation factor. When this technique was reinstated in the late 
I980’s using the intensified array detector and a modified electronic lamp power supply, the 
radiative transfer model required a factor, I, = 22, (relative to coded values in Fig. C-2) to 
accurately match laboratory- measured transmittances down to 0.2. Without ! t . the model 
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predicted transmittances faLrly accurately down to about 0.8 but underpredicted at lower 
transmittances with increasing deviation. Since I t had little effect on transmittance calculations 
above 0.8, a' was determined from a wide range of laboratory conditions (P and T) 
approximately independent of I t . The results for a' were consistent with the 1981 reported 
value. Also, a two parameter non-linear least squares fit to all the laboratory data with the 
temperature exponent, n, fixed at 1 .5 gave L » 22 and again a' consistent with the 1981 values. 
The study concluded that radiative transfer model calculations and laboratory data were in 
great agreement over the full range of laboratory conditions for I c = 22 and a ~ 3I800*P/1 15 . 
In recent years, the broadening parameters for individual lines have been measured directly 
(Ref. C-7) and reported as a' = 6927*P/T U5 . It is both important and gratifying that this agrees 
extremely well with a' = 3180Q*P/T M . Using the Ref. C-7 value in the code, I c - 22 was still 
the best fit value. Relative to the Fig. C-2 lamp characterization values, I c and a' were about 
equivalent for both the (0,0) and (0,1) bands. 

The instrument-dependent code calibration parameter, a', was determined independently 
for the continuum absorption system using continuum lamp transmittance measurements over a 
similar range of laboratory measurements. For the continuum absorption system, the broadening 
constant, C * 23250, gave the best overall fit to laboratory data with the temperature-dependent 
exponent, n, fixed at 1.5. 

A radial inversion scheme, referred to as an onion peel, was used to determine a radial 
profile of NO number density from a scries of LOS measurements made at several radial 
locations, as illustrated in Fig. C-3 during a steady state-test condition. The flow field was 
modeled as concentric “homogeneous” zones, one for each measurement. The density in the 
outer zone was determined using the outermost measurement. Stepping inwardly, a new zone NO 
density was determined for each successive measurement. The static temperature and static 
pressure along the optical path were inferred from probe measurements as discussed in Appendix 
A. Probe measurements at similar radial positions on opposite sides of the plume were averaged 
and radial symmetry then assumed. Measurements of ambient test cell temperature were assigned 
to positions beyond the lip of the diffuser, and ambient pressure assigned to radial positions just 
beyond the nozzle exit radius, Re = 30 cm. Radial distance was normalized to the nozzle exit 
radius so that the extent of the engine exhaust was; R/Re as 1.0 or a little beyond, depending on 
exhaust expansion for each test condition. The probe measurement grid was sparse, and thu,s 
interpolation was required to determine the path average temperature, pressure, and axial flow 
velocities for each homogeneous zone along each LOS measurement path. Averages were area 
weighted over the width of the optical beam within respective zones. Figure C- 4 shows 
representative plots of zone-averaged temperatures, pressures, and axial velocity for an LOS 
through the center of the exhaust. Probe-inferred values used io determine the zone averages are 
superimposed on the plot. These data were acquired for the engine test condition described by 


81 


AEDC-TR-96-3 


the combustor inlet temperature T3 = 733 K at 9.1 -km altitude. Refer to Table 3 in the main text 
for summarized engine/facility test conditions. The zone path segments are obvious in the plot 
with approximately 2.54-cm segment lengths from the center to just outside the exhaust flow 
field, and about 5.08-cm segment lengths in the test cell gas-recirculation region. Zones beyond 
the nozzle exit radius represent measurements through the test cel! air, which contained small 
quantities of recirculated exhaust gases. 

€.4.0 TEST CELL INSTALLATION AND MEASUREMENTS 

The relative placement of the NO-UV instrumentation and optical beam paths relative to 
other measurement and facility components in the test cell are illustrated in Fig. 2 of the main 
body of the report. Preliminary measurements during engine checkout verified that recirculated 
exhaust gases containing NO would be present at low levels in the test cell. Since even low 
levels of NO over rather long paths can have a significant effect on the results, measurements 
were also performed outside the exhaust flow field to quantify recirculated NO concentrations. 
Quantification, pf course, assumed cylindrical symmetry for recirculated gas components 
surrounding the exhaust flow field. 

The source radiation was transmitted through a 9-m length, 400-pm-diam, fused silica 
fiber to a collimating lens mounted onto a traversing table near the engine exhaust. The radiation 
was imaged through the exhaust to an aluminum-surfaced mirror (mounted on a second 
traversing table) and reflected back through the exhaust to a second lens-fiber combination 
mounted on the traversing table. The beam diameter was approximately 2 cm. The transmitted 
and reflected beams were within the first 12 cm of exhaust from the nozzle exit plane. The 
resonance lamp beam was about 6.35 cm above that of the continuum lamp. The reflected 
radiation was transmitted to (he entrance slit of the spectral detector system. All optical lenses 
and fibers were made of UV grade fused silica. 

The slide tables traversed the collimated beam across the flow' field in a plane parallel to 
the nozzle exit plane, as illustrated in Fig. 2 of the main text. Each set of traversed data, for the 
resonance absorption system, consisted of measurements through the purged path at the top of 
the diffuser (with and without the shutter closed) and measurements at fixed radial positions 
from R '= 70 to 36 cm in 5.08-cm increments and from R - 36 to 0 cm in 2.54-cm increments. 
The set of continuum absorption measurement positions were offset about 6.35 cm below the 
resonance absorption measurements, except, of course, for measurements through the purged 
tube. A traversing system output voltage was monitored by each OMA system for posttest 
determination of memories (spectra) acquired while the tables were in motion 

During an earlier engine checkout run, measurements with the lamp shuttered were found 
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to be indistinguishable from background signal levels. Therefore, acquisition time was 
significantly reduced by not shuttering the lamp except at the centerline position. These 
centerline data were used to verify that exhaust radiation levels were negligible for each engine 
condition. It is also significant that the measurements were not affected by ambient lights within 
the test ceil. Background levels were primarily integrated detector noise as discussed in OMA 
operation manuals. 


C-5.0 DATA REDUCTION 


Spectral lamp transmission measurements, I T (X,), at a particular LOS (radial) position were 
reduced to spectral transmittance according to the equation 

iA>-W MV- W /ro) 

l«(X ; )-\(X,)”l n (\)-l b {^) 

where X, denotes the spectral center of the \ ih element of the detector array. loiX,) and J h (X f ), 
respectively, are measurements of lamp reference and background at the purged path location. As 
mentioned earlier, measurement of lamp-shuttered signal levels at the exhaust, ^(X.,), were 
indistinguishable from background, so I,.(X) can be approximated by 1 S (X,). Lamp and 
transmittance spectra (background subtracted) are presented in Figs. C-5 and C-6 for the 
resonance and continuum absorption systems, respectively. Prominent NO gamma band regions 
are labeled in the figures. Absorption features at the (0,0) band are readily evident, but further 
corrections were necessary to remove position-dependent lamp signal variations and attenuation 
effects of the exhaust: and gas recirculation. These effects are exhibited as overall attenuation at 
spectral regions away from the (0,0) band absorption feature. A curve fit through spectral regions 
just above and below the (0,0) band feature served as the base line to quantify (0,0) absorption, 
as illustrated by the “corrected” curve in Figs. C-5b and C-6b. Thus, the transmittance at the 
second bandhead of the NO gamma (0,0) band can be read directly from the "corrected” curves. 

Formerly, spectral bandhead transmittanees were determined using pretest lamp reference 
spectra acquired at each radial position and then corrected for test media attenuation effects. Data 
were reduced both ways and gave equivalent results; i.c., differences in second bandhead 
transmittances from the two processes were far less than statistical uncertainties. The purged- 
path lamp spectra allowed an independent measure of lamp reference per data set and all the 
advantages that go with measuring lamp reference spectra within seconds of acquiring 
transmission measurements. Statistical uncertainty is discussed Inter. 

Resonance absorption system radial profiles of (0,0) second bandhead transmittance, 
acquired at the five engine power settings at 9.1 -km altitude, arc presented in Fig. C-7a. These 
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profiles were smoothed using a linear 3-point smoothing algorithm before insertion into the 
radial inversion scheme. The two bottom curves, almost coincident, show smoothed and 
unsmoothed transmittance profiles for the same engine test point, illustrating representative 
affects of the smoothing algorithm on these profiles. Although random-like variations were 
reduced, changes in transmittance values were generally much smaller than statistical 
measurement uncertainty. Multiple sets oi transmittance profiles were acquired at engine 
conditions held at steady state for several minutes. Several sets of transmittance profiles acquired 
at the same, engine condition were extremely repeatable with variations less than statistical 
measurement uncertainty. Multiple transmittance profiles acquired during a steady-state test 
condition were averaged before application of the radial inversion. 

Continuum absorption system radial p ottles of second bandhead (0,0) transmittance, also 
acquired at the five engine power settings at 9.1 -km altitude, arc presented in Fig. CT7b. 
Comparing these profiles with resonance absorption profiles in Fig. C-7a, it is obvious that the 
continuum absorption system was less sensitive to equivalent NO density at the same test 
conditions. Additionally, continuum absorption was not observed for several low-power engine 
settings which produced low NO densities Preliminary results on selected data sets indicate 
higher uncertainties than results obtained from the resonance absorption data. The higher quality 
and more complete resonance absorption data were analyzed and reported herein. The continuum 
absorption data will be further reduced and analyzed for publication in a later report. 

C.6.0 ANALYSIS AND RESULTS 

The radiative transfer model, Eq. (C-l), incorporated into the onion -pee I radial inversion 
scheme was applied to radial transmittance profiles for each steady-state engine test condition. 
Static temperature and pressure profiles, required for the radiative transfer model, were inferred 
from probe and facility measurements and averaged along “homogenous” zones for each LOS 
optical path. T his measurement technique determined NO number density (molecules per 
volume) profiles. For a more direct comparison to probe sampling measurement profiles, radial 
NO-UV number density (cm *) profiles were converted to volumetric fractions (parts per million 
by volume, ppmv) using the ideal gas law. Profile comparisons are presented in Figs. C-8 
through C-10 for the same combustor temperature setting, T3 ~ 7 J3 K, at three different 
altitudes. 3.! (SLS), 9.1, and 12.2 km, respectively, The NO-UV densities tend to be similar in 
shape to the sampling, but lower by 10 percent for 9.1 and 12.2 km, and about 18 percent for 
SLS. The stronger deviation might be expected at higher engine power settings for the SLS 
condition because of the combined effects of greater quantities of NO recirculation in the test cell 
and possible non-radiul symmetries (especially top to centerline) in the recirculation region. 
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The more direct comparisons of density profiles were limited to the few test conditions for 
which spatial probe sampling profiles were measured. For comparisons over the full set of test 
conditions, emission indices were calculated from NO-UV resonance absorption density profiles, 
aerodynamic flow-field properties derived from the rake data, and fuel flow rates measured by the 
test facility. An NO emission index EI(NO), using the molecular weight of NO : , was determined 
for each engine condition according to 


Fil(NO) = 


1000* M NOj 
N A *( fuel ilow rate) 


2£( n No)> v X]Aj 


(C-3) 


where N A is Avogadro's number and M no ^ is the gram molecular weight of NO,. The quantities 
n N0) and VXj are average zone NO number density and exhaust velocity, respectively, and is 
the zone area. These Ei(NO) values are presented here in Table 5 and in Figs. 13-17 in Section 
4.0 of the main body of the report for comparisons to probe sampling and infrared tunable-diode 
laser technique results. 

UV EI(NO) values were in good agreement with probe sampling data, matching well 
within respective uncertainty bands. However, NO-UV EI(NO) values were usually lower 
throughout the test. This is consistent with the density profiles compared earlier. A contribution 
of uncertainty not yet addressed in these comparisons relate to aerodynamic flow-field properties 
determined from rake data. As mentioned, total temperature and Mach/flow angularity probes 
wore sparsely spaced giving a coarse definition of property profiles. The slow varying nature of 
these properties across the exhaust flow field gives good confidence in values reported, except 
near the edge where changes are more abrupt, especially in velocity. The product of velocity and 
large zone area, near the edge can result in a variation of a few percent of EI(NO) in Eq. (C-3) 
calculations, for just centimeters uncertainty in the radial position for which velocity drops 
abruptly. 

An engineering s atistical analysis process, described in Ref. C-8, was applied to steady- 
state test conditions for which several repeat measurements were performed to give confidence 
bands for the NO-UV resonance absoiption results. Lamp reference and position-wise lamp 
transmission spectra were treated as independent measurements with random statistical 
uncertainty. The variance-covariancc matrix developed for a set of profile transmittances took 
into account codependency on lamp reference spectra, curve-fit correction factors, and the linear 
smoothing algorithm. Measurement uncertainty and a matrix of partial derivatives of NO number 
density with respect to transmittance for each zone of each LOS measurement path produced one 
sigma variances for NO number densities at each radial location per set of data. 
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Independently, statistical analyses were performed on the multiple data sets acquired at 
steady-state engine condition.'-, treating each set of spatial transmittances as independent 
measurements with random uncertainty. The resulting NO number density uncertainty was 
similar, tending to validate assumptions used for the lengthy and tedious engineering statistical 
analysis process required for single data sets for many of the engine conditions. 

The position-wise number density uncertainties within the exhaust flow for several test 
conditions investigated ranged from ± (6 to 18) percent, with the higher uncertainties relating to 
data at upper T3 values for the simulated sea-lcvel-static (SLS) condition. Statistical iincert;i ,n ty 
on position- wise densities determined for the recirculation region were on the order of i 15 to 
percent. Statistical uncertainty for emission indices r mgsd loan ± (7 to lb) percent assuming 
accurate velocities from the rake data. 


Fo- completeness, resonance absorption was also observed at the (0,1) band for a couple of 
the higher temperature conditions. Absorption at this band and the (0,0) band will allow 
simultaneous determination of NO density and ground-state vibrational temperature (which in 
this case should be equivalent to kinetic gas temperature) in the inversion process. Results of this 
analysis will be detailed in a later report. 

C-7.0 SUMMARY AND CONCLUSIONS 

NO-UV resonance absorption measurements were performed at all test conditions 
throughout the test period consisting of several engine power settings at simulated altitudes 
ranging from SLS to 15.2 km, Spectral transmission measurements at multiple radial stations 
across tiie exhaust flow field and surrounding region containing low-level recirculated exhaust 
gases were radially inverted to provide radial profiles of NO number density without prior 
assumptions on radial species concentration profiles, required for analysis of single !ine-of-sight 
optical measurements. Radial NO number density fenr 1 ) profiles were converted to radial 
profiles of concentrations (ppmvj, and subsequently to emission indices, using radial profiles of 
static temperature, static pressure, and exhaust flow velocity inferred from probe rake data. The 
NO emission index is a single value for an engine test condition. 

Tiie NO-UV measurement system proved fe he reliable and robust. Measurements were 
extremely repeatable during steady-state test conditions, indicating temporal stability of engine 
and test facility parameters. 

NO-UV radial concentration profiles were compared with radial profile gas sampling data 
and agree well (within unvCrtainty limits) with probe sampling measurements, for both the 
exhaust flow field and the surrounding gas-recirculation regions. The NO-UV concentrations 
tended to he lower than probe sampled levels within the exhaust region. NO emission indices 
determined from NO-UV profile concentrations deviated only a few percent from the probe 
sampling results. 
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The NO-UV resonance absorption technique performed well throughout the test program, 
produced results consistent and in good agreement with probe sampling and infrared tunable- 
diode laser techniques, and is recommended for future test measurement programs. 
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Figure C-7. Radial profiles of (0,0) transmittance a»t 9.1-km altitude. 
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APPENDIX D 

INFRARED TUNABLE-DIODii LASER SYSTEM 

Joda Wormhoudt 
Richard Miake-Lye 

D-l.® TDL INSTRUMENT 

Tunable lead salt diode lasers have a long history of spectroscopic applications (Ref. D-l), 
and have recently been applied with considerable success to sensitive detection of atmospheric trace 
gases (Ref. D-2). Our system had been applied to combustion gas flows twice before, first in a 
combustor simulator at NASA Lewis Research Center (LeRC) (Ref. D-3) and then in a short engine 
test at AEiDC (Ref. D-4). 

Our diodes are housed in a liquid nitrogen dewar, with their nominal operating temperature 
controlled by a small resistance heater. Coarse frequency tuning is achieved by changing the diode 
temperature and thereby its refractive index and cavity mode frequencies. Fine tuning is done by 
scanning the injection current, which also has the effect of slightly changing the diode temperature. 
This current scanning can be done rapidly, allowing the averaging of many spectra whose spectral 
features tire unperturbed by low-frequency fluctuations in the transmission. In the work reported 
here, the scanning rate and number of spec tral points acquired was such that the scan of each diode 
took 0.5 msec. 

In many applications, simultaneous measurement of concentrations of two molecules 
provides much greater understanding than measurement of either one, or even sequential measure- 
ments of both. Diode lasers emit light in a restricted range of frequencies, so that different laser 
diodes are usually required to delect each molecule of interest. Our typical cryogenic dewar 
configuration houses four laser diodes behind an infrared (IR) transmitting window, with 
translation of the dewar se rving to select which diode emits into the collecting optics and so through 
the entire optical train, until laser light is focused onto a detector located in the same dewar. We 
therefore dre w the conclusion that the usefulness of the TDL instrument would be greatly increased 
if data could be acquired using two diodes at the same time, and set to work reconfiguring the op- 
tical layout of the laser table into the arrangement shown in Fig. D-l . 

Our focus on in- situ measurements in aircraft engine exhausts motivated two other choices 
in this design. First, we could gain space on the optical table by eliminating the long path absorption 
cel) and therefore the extractive sampling capability used in the NASA LeRC field test. Second, our 
measurement time scales were such that it was adequate to use a single detector, with discrimination 
between the two lasers achieved by alternating lascr-on (frequency scanning) and laser-off periods. 
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Therefore, the optica) layout drawn in Fig. D-l is capable of collecting light from two laser diodes 
(one in each of two- banks of four on opposite sides of the dewar), sending two coincident, 
collimated beams to a retro re fleeter, then collecting the return beams anci focusing them onto the 
same detector element. 

We begin a discussion of the details of the optical layout by tracing the paths of the main or 
diagnostic beams. Infrared light from the laser diodes is collected by reflective microscope 
objectives (15x) and focused onto 200-jim pinholes, which define input apertures. The pinholes are 
used only during alignment, so they are mounted on removable kinematically indexed bases. The 
microscope objectives are mounted on X-Y-Z translators to allow accurate focusing into the fixed 
apertures. Past the input aperture, curved mirrors create approximately collimated beams. The two 
laser beams are then combined at the main beamsplitter, and the co-aligned beams pass to a pair ot 
flat mirrors and onto a large gimbal-mounted steering mirror which directs the beam off the table 
and into the engine test cell. The return beam is collected by the large steering mirror and a curved 
mirror of similar diameter, passed through several more flat mirrors and finally focused onto the 
main beam detector. 

The remaining combined beam associated with the main beamsplitter goes into the reference 
optical path, where it may be sent through a gas calibration cell and/or passed through a grating 
monochromator, then focused onto the reference detector. The calibration cel! contains a sufficient 
pressure of the gases of interest (in this case, NO, NO v and N 2 0) that their strong absorption 
features fully absorb the laser light. This makes it easy to recognize that the laser is operating in the 
right; wavelength region, and to check for the presence of other laser modes. These strong 
absorption features can also be used to set the frequency scale of the laser scan. The monochromator 
is used to provide a coarse measurement of the laser frequency during setup of the instrument, and 
to determine whether the diode is emitting at different wavelengths (different modes). Since no 
monochromator is available in the main beam path, diode operating conditions must be chosen that 
allow single mode operation. 

There is a parallel visible optical system for alignment and setup. A red HeNc ‘‘trace laser 
beam passes through the dichroic beamsplitter downstream of each microscope objective, and is 
coaligned with each infrared beam. Coalignment is guaranteed by focusing the beam through the 
input aperture. The trace beam is an indispensable aid for alignment of the optical system. In 
addition, the trace beam is used for accurate calibration of the monochromator via higher order 
diffraction. The fourth port of the dichroic beamsplitter can be used to observe the laser diode 
during alignment. An eyepiece mounted at the position conjugate to the pinhole forms an effective 
150x microscope. 
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We set up the instrument with an NO diode on one side of the dewar and N0 2 and C0 2 diodes 
on the other, with the C0 2 diode operating at the beginning of the test. At several points during the 
test, the collection optics were moved to pick up light from the N0 3 laser, then moved back to the 
C0 2 laser. The change from one diode to another followed by optimization of transmitted laser 
intensity could be successfully completed in a few minutes. 

The temperature and current of each laser are controlled by a Laser Photonics controller. All 
inputs to this controller were made via a program running on the data acquisition aH analysis 
computer. In a typical experiment, the laser temperature is held constant while the current through 
the laser is modulated with a computer-generated sawtooth to sweep the output frequency across 
the IR transition- The sawtooth was generated in 150 discrete steps for each laser by a digital-to- 
analog converter board. The total of 300 display points in the two laser sj>ectra is swept at a rate of 
1 kHz. During approximately 10 percent of the duty cycle, the laser current is dropped below the 
lasing threshold to provide a precise measurement of zero light intensity. 

The data acquisition techniques., involving a Scientific Solutions Lab Master analog-to 
digital converter hoard and direct memory access transfer to the extended memory of the data- 
acquisition computer, remain as described previously (Ref. D-3). Although the TDL instrument 
could in principle produce real-time records of gas concentrations, this would require input of 
probe rake data (or assumptions about properties along the line of sight (LOS) from other 
sources). Our data acquisition therefore consisted of the recording of a large number of spectral 
files covering segments of both the NO and CO : (or N0 2 ) regions. These were recorded as 
averages over periods of 10 to 60 sec. 

D-2.Q DIAGNOSTIC I MPLEM ENTATION IN THE ENGINE TEST CELL 

Figure D-2 shows two views of how the TDL diagnostic was applied to the engine exhaust in 
the AHDC altitude test cell. The TDL instrument was placed on a table directly below the optical 
access port into the test cell. The TDL beams were directed into the test cell by a 45-deg mirror 
mounted on the outside of the access poll flange. A 0.95-cm thick CaF 2 window with a I -deg wedge 
(to suppress interference fringes) was mounted on the access port using a retaining plate with a 10- 
cm-diam clear area. Because the access port was offset from the engine exit plane location, the 
beam next: passed th rough a 1 50-cm long aluminum tube containing two flat mirrors, This peri scope 
was purged with gaseous nitrogen to protect against both recirculating exhaust, and oil and water 
mists from the engine, and the exhaust treatment device called a diffuser. The TDL beams wore 
returned, horizontally displaced by about 2 cm by a retrore fleet or. The HcNc trace beam could be 
followed through the entire optical path, which allowed precise aiming of each steering minor, 
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The potential dangers to our mirrors of the diffuser water spray and engine oil lumes 
mentioned above concerned us enough that an arrangement was made to place baffles in front of 
the periscope and retrorefiector when the TDL diagnostic was not operating. This was done by 
attaching aluminum baffle plates on the AEDC ultraviolet cross-flow absorption device, k 
addition, purged tubes extended from the periscope (50 cm in length) and from the retroref.ector 
{45 cm) toward the exhaust flow, leaving an unpurged length along the LOS of about 238 cm. 

B-3,0 EXAMPLES OF OBSERVED SPEC TRA 

Figure D-3 presents two examples of raw diode intensity spectra for the TDL diagnostic 
beam, one taken while the laser beam was passing through the exhaust and the other while the beam 
path was blocked by the UV instrument. The vertical axis is simply intensity at the detector 
expressed in mV, with (he lowest point on the beam-blocked spectrum being shifted to zero, and 
the exhaust spectrum also shifted until the laser-off points at the right-hand side of the figure 
coincide with the beam-blocked values. The horizontal axis is in display units or data acquisition 
points over the laser scan. The NO absorption line used in almost all measurements is clearly seen, 
with the wing of the first of two neighboring water absorption lines contributing to the decrease in 
transmission at the right-hand side of the figure. These water lines have fairly high lower-state 
energies and hence are substantially stronger in the hot exhaust plume than they are at ambient 
temperature, even over the much longer optical path outside the test ceil. 

The NO, scan presents smaller absorption features than those seen in the NO scan, but it will 
be seen below that the absorption lines indicated are due to known lines of NO, and H,0. A mode 
change is seen at about 20 display units. Observation of the reference leg intensity as the 
monochromator grating is scanned confirms that the diode laser intensity below this point was in 
an entirely different wavelength range than in the rest of the scan. The diode is emitting in only one 
mode throughout the entire scan, but there is an abrupt change from one mode to another. (Clearly, 
in insidious cases there could be no telltale discontinuity, and only the monochromator diagnostic 
would make the situation clear.) In Fig. D-4, wc present a corresponding raw intensity spectrum at 
the reference leg detector, while the. reference laser beams are passing through a short cell contain- 
ing substantial amounts of NO and NO... In the NO scan, two strong NO line features arc seen. The 
line at 222 display units is the diagnostic NO line indicated by the arrow in Fig. D-3. The second 
line is exactly coincident with the lower frequency of the two water lines in Fig. D-3. Therefore, it 
is of no interest for the measurement of NO in exhausts, but it would have to be taken into account 
if one wanted to use the water lines for high accuracy measurement of the water mixing fraction in 
the plume. 

Turning to the NO, laser reference scan, we see the high density of NO- lines which adds to 
the difficulty of sensitive detection of this molecule. Furthermore, it will be seen that a similar den- 
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sicy of water lines also exists in this region. The line we will use as the diagnostic line is the strong 
line near the center of the N0 3 scan, at 77 display units. 

D-4.0 DATA ANALYSIS PROCEDURES AND ERROR LIMIT ESTIMATES 
D-4.1 DATA ANALYSIS 

Several steps are required to convert the raw diode intensity scans into transmission spectra 
and deduce molecular mixing fractions, mass fluxes, and emission indices. In the first step, 
absorbances of the diagnostic iines are quantified by least squares fits, using the TDL data 
acquisition program in data analysis mode. This is the routine that allows real-time fitting and 
concentration measurement, here applied to stored spectra. Lines are fit to Voigt line profiles using 
a nonlinear least squares rout ine based on the Levenberg-Marquardt approach (Ref. D-5). The base- 
line, varying from unit transmission due to diode laser power fluctuations and residual etalon 
fringes, is represented as a slowly varying polynomial of low order. Lines or sets of lines in six 
speck'd segions were fit: the NO line pointed out in Fig. D-3 (near 1897 enr'), the pair of water 
lines to its right, an NO line near 1909 cm -1 (for two engine conditions), a pair of N0 2 lines, one of 
which has an overlapped water line (ail near 1597 cm -1 , observed during six engine conditions), and 
a CO, line near 2225.5 cirr 1 (observed during one engine condition). 

Figures D-5 and D-6 show two of the spectral regions that were analyzed. In each figure the 
solid line is the observed trace as converted to a transmission scale by the TDL data analysis pro- 
gram, through its fit of a polynomial baseline (this calculation also uses as an input the relative po- 
sition of the laser-blocked intensity trace shown as the dotted line in Fig. D-3). The ratios of 
Doppler and pressure-broadened line widths in the fitting models were constrained, using HfTRAN 
parameters (Ref. D-6) and a single representative temperature and pressure that changed with en- 
gine condition. In the fitting of the 1897-cnrr : H 2 0 doublet, there is an additional constraint of a 
fixed intensity ratio and an additional fitting parameter, the spacing between the two lines. In the 
N0 2 region, all lines contributing to the two features in Fig. D-6 are included, with NQ 2 and H 2 () 
concentrations varied separately to achieve the best fit. 

The second step is to convert the optical depths (absorption line peak intensities) determined 
by the least squares fits into a set of column densities over the nonuniform optical path through the 
exhaust plume. Our approach uses predicted spectra generated from the HJTRAN database using a 
program written at Aerodyne Research, Inc. (AR1), which takes sets of temperatures, pressures, 
species mole fractions, and path lengths for each segment of the nonuniform LOS through the 
exhaust plume and calculates a total absorption spectrum. The result of this model calculation is a 
line peak absorbance to be matched to the observed value from the least squares fitting discussed 
above. The sets of temperatures, pressures, and path lengths used in our analyses were those used 
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in the analysis of the AEDC NO resonance lamp data. This model assumed the plume to be made 
up of 20 concentric annular cylinders with constant properties equal to linear interpolations of the 
probe values. The final set of parameters we require, the relative mixing fractions in tit v e 20 
cylinders, can be obtained from several sources, and we used a number of sets to evaluate our 
sensitivity to the assumption of a particular model mixing profile. Examples of two such mixing 
models are shown in Figs. D-7 and D-S for two engine conditions. 

One mixing fraction profile is derived by taking ratios of differences of probe total 
temperature measurements and an ambient total temperature (converted from static temperatures 
outside the exhaust flow, or taken from probe total temperature values outside the exhaust) to that 
difference at the plume center. This relative mixing fraction is near unity for the innermost cylinder, 
and decreases fortheoutercylindcrs. The other mixing fraction show n in the figures is derived from 
the inversion of the AEDC NO resonance lamp data sets, and consists of NO mole fractions 
normalized to unity at the centerline. The absolute value of the centerline mixing fraction is then 
determined as the multiplicative factor that brings the model calculation of line peak transmission 
into agreement with the value determined from the least squares fit of the ohserved line. The 
resulting model spectra for the two models described above are averaged to give the dotted lines in 
Fi» >. D-5 and D 6. Figure D 6 also contains mode! spectra for the individual contributions of NO* 
and H 2 0. In order to plot these model calculations on the same scale as the data and TDL analysis- 
program fits, the line positions and frequency scales were set by the reference line positions shown 
in Fig. FM. 

D-4.2 DERIVATION OF CONCENTR ATION AND MASS FLUX VALUES 

The sampling probe not only produced spatially resolved values of temperature and pressure 
but also velocity. These values can be used to convert our single LOS measurement to a total species 
flux in the exhaust. Our procedure was to use the sa n\c 20-cylindcr model, multiplying cylinder 
areas, species number densities (fanned from mixing fractions and total densities computed from 
static temperatures and pressures), and velocities with conversions involving the molar volume and 
the molecular weight to yield a mass flux in g/s. In the aircraft emissions community, the results of 
such mass flux measurements are typically expressed as emission indices (HI), the ratios of g/s of 
the exhaust species to fuel flow' rate in kg/s. Our conversion from mass flux tu HI was made using 
the tabulated fuel flow rates. 

As a check on this procedure, we performed the same mass flux calculations except with the 
mixing fraction omitted. This total mass flux should match the tabulated input mass flux, a quantity 
measured at a level equivalent to the fuel flow rate mentioned above. It turns out that in most (but 
nm: all) cases, the total mass flux computed from probe rake measurements is substantially less than 
the inlet mass flux. Discrepancies are not surprising, given (he coarse grid on which probe rake 
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sampling measurements were made. However, it is not possible to be certain whether a correction 
factor that converts probe rake integrations into inlet total mass fluxes should be applied to exhaust 
species emission indices as well, If the integration of probe rake values has underestimated 
contributions to the total mass flux from regions near the outer boundary of the exhaust flow, the 
corrections to an NO or CO ; emission index will be much smaller than the correction to total mass 
flux, because the mixing fraction of exhaust species is going to zero in this region. On the other 
hand, corrections stemming from inaccuracies in flow parameters near the centerline will apply 
equally to the total mass flux and to exhaust species emission indices. Because of this uncertainty, 
we did not apply individual corrections factors, but used their average to estimate the magnitude of 
errors involved in using probe rake data to convert concentration profiles to emission indices. 

The other contributions to an error estimate which we considered were the statistical errors 
in sets of multiple optical depth measurements, and an estimate of the systematic errors involved in 
conversion of optical depths to concentration profiles. The former actually involve two classes of 
nominally comparable measurements. The first class involves multiple measurements while the 
engine remained at a given operating condition. Analyzed in the same way, multiple spectra 
typically gave remarkable repeatable optical depths, with standard deviations of a few percent 
except for low signal idle conditions. The second class of multiple measurements involves variation 
of the order of the baseline polynomial used as part: of the TDL data analysis program least squares 
fitting to determine line peak optical depths. Polynomials cf order three through five were used. 
Again, differences were typically small except in low signal cases. In some of those cases, there 
may be interfering lines or etalon fringes, arid additional analysis of individual spectra may reduce 
these error limits. 

Our estimate of the errors involved in converting optical depths to concentration profiles (that 
is, assigning a peak concentration scaling factor to an assumed mixing profile within a model 
calculation of the absorption spectrum along a nonuniform path) focused on the differences 
resulting from different mixing profiles. (Example calculations performed by varying the assumed 
temperature profiles, for example, indicated lower sensitivities to the other model parameters.) The 
model calculation discussed above, of predicted optical depth for a given set of model parameter:; 
including mixing fraction and of conversion factor (independent of spectral region) between 
centerline mixing fraction and emission index, can be combined into a single conversion factor, 
which multiplies peak optical depth for a given spectral line to convert it to emission indices. In 
Table D- 1, these conversion factors are listed for the 24 steady engine conditions and the NO 
absorption line near 1897 cm for the two mixing model examples shown in Figs. D-7 and D-8. 

Tbc.se two models, while broadly similar, differ in at least two significant ways. First, while 
(he mixing profiles based on ratios of total temperature differences almost always have their peak 
values on the centerline, the IJV NO profiles typically show a peak in mixing fraction at an 
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intermediate radius. Although uncertainties in the inversion of UV lamp data increase for the lines 
of sight near the centerline, wc believe that the indication of a dip in the middle of the concentration 
profile is qualitatively correct, since for the three cases where spatially resolved sampling probe 
measurements were made* such a central dip ir clearly and consistently observed. The significance 
of such differences in profiles is that an outer cylinder makes the same contribution to the observed 
line-of-sight optical depth as one near the centerline with the same column density (product of 
number density and path length), but its contribution to the mass flux and hence the emission index 
is larger by the ratio of the radii. 

The second difference, of significance in the sea-level cases such as the one shown in Fig. 
D-7, is that our particular implementation of a model based on ratio temperature differences 
accounts for very little recirculating exhaust gas outside the jet radius, while the UV NO based 
profile in Fig. D-7 indicates a significant concentration. However* during this case and a few 
others* sampling measurements were made at three points outside the jet radius (thtr outermost at 
just inside 3 radii). These measurements indicate that the true exhaust mixing fraction in this 
region along the centerline LOS lies between the values indicated in Fig. D-7. It must be kept in 
mind that the IJV NO lines of sight that determine the values shown in the dotted curve in Fig. 
D-7 are above the plume, not on either side* and it may be that more recirculating gas is found 
near the top of the test cell. Also* a third full set of models was used to analyze the data* with the 
mixing fraction determined by ratios of static temperature differences. This is clearly not 
defensible on first principles, and indeed gives a mixing fraction within the exhaust flow that is 
much too peaked at the centerline compared with any other method of estimation. However, this 
model has the advantage that the recirculation sampling measurements are ail fairly well 
reproduced (most to within a factor of two), so it can function as a correlation of expected 
concentrations in the recirculation region for all cases. Inclusion of these conversion factors in 
our evaluation did not increase the estimate of error limits due to systematic errors in the assumed 
mixing profiles. 

Numerical values resulting fiom the above three estimates of contributions to an overall error 
limit arc listed in Table D-l . In that table, NO emission indices (from 1 897-cnrr 1 line observations) 
and the resulting estimate of a total error limit, in El units* are given in the second and third 
columns. The fourth column, the fractional standard deviation i.n multiple measurements of 
comparable optical depths* is seen to vary from less than 1 percent vin one idle case) 'o more than 
30 percent. The next two columns arc the absolute values of the calibration factors which convert 
peak optical depth to EI(NO) within the two models discussed above. The deviations from unity of 
the ratios of these two columns are averaged to give an estimate of this uncertainty level as a 
fraction of Id. This average is then applied uniformly to all cases. The same procedure is applied to 
the last column, which is the ratio of total mass flux measured at the engine inlet to mass flux 
deduced from an integration of exhaust probe rake values. The fact that this converts a correction 
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factor, which in almost all cases increases the reported emission index into an error estimate which 
implies positive and negative corrections are equally likely, is not a problem, in our opinion. There 
are other potential errors in our model-based analysis which could result in systematic downward 
corrections, so our choice of expanding the error bars rather than employing the correction factors 
seems reasonable. When the error estimates derived from averages of differences from unity of 
mixing model conversion factor ratios and total mass flux ratios are added in quadrature, the total 
fractional error limit assigned to these systematic errors is 19 percent. This is added (in quadrature) 
to the statistical error estimates in optical depth listed in the fourth column to give a total fractional 
error, which is then multiplied times the EI(NO) to give error limits in El units, tabulated in the third 
column, and used in the comparison plots in the main body of this report. 

Our analyses of the smaller number of data points obtained in other spectral regions indicated 
that similar 20- to 30-perceni error limits in emission indices should apply. One final comment can 
be made based on the two columns of HI conversion factors in Table D-l . If the idle condition points 
are excluded, it is interesting to note that all other conversion factors vary' by not much more than 
a factor of two, even though conditions vary over a wide range of power settings and simulated al- 
titudes from sea level to 15,2 km. In other words, although mud) of our discussion has focused on 
how the cm rent TDL technique is dependent on characterization of the optical path for the 
derivation of accurate values, this shows that factor-of-two emission indices could he obtained 
directly from optical depths with only minimal corollary information about the exhaust plume, such 
as previous observations or a CFD calculation. In the present test where the focus was on 
comparison to state-of-the-art NO measurement techniques, such approximate determinations are 
of no interest. However, it is possible to imagine situations in which the TDL technique could be 
extended to new species, for which even order-of-magnitude measurements arc not now available, 
and used in a screening mode to identify which species, engines, and operating conditions require 
more careful and quantitative studies. In those situations, factor-of-two values without the need for 
probe rake instrumentation might be very useful. Of course, we hope to develop TDL techniques 
which use line pairs (such as an NO and a C0 2 line) to make self-calibrating measurements, as is 
currently done by deriving sampling Ei(NO) measurements from the ratio of NO and C0 2 
concentrations. However, some additional development work will be required before wc can be 
sure this can be done, even for the species wc already detect. 

The full set of observed emission indices obtained to date is given in Table D-2. The 
generally good agreement between the KI(NO) values and those from the two other measurement 
techniques has been discussed in the main body of the report, where comparison plots showed 
identical trends in almost all cases and agreement within the estimated error limits, again in almost 
all cases. A few anomalous points can be mentioned: the 9.1 -km idle NO observation is clearly too 
low, whereas at 15.2 km, differences at one or more power settings lend to the appearance of a 
different trend than observed by the other techniques. It may be that additional analysis of the 
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spectra in these cases will resolve some of these anomalies. For two engine conditions, a second 
line at 1909 enr 1 was used to measure NO in addition to the line at 1897 cm 1 . The NO 
concentrations determined from the two lines were in good agreement, as expected. The N0 2 and 
H 2 0 values obtained in the 1 597-cnr : region are dearly sensitive to assumptions about the baseline 
position (in this spectral region, lines are sufficiently overlapped that the baseline is not directly 
observed). Still, with the exception of the 9.2-km idle point, the water emission indices agree with 
the theoretical value of 1233 to well within the expected uncertainty limits. The water lines in the 
1897-cnr 1 region arc inadequate for accurate measurements in several respects, so we evaluated 
one observation, found {not surprisingly) a substantial disagreement with the theoretical value, and 
do not plan further investigation of this region. The C0 2 observation is the only one in which all 
lines in the scan are in the HITRAN92 tabulation. Spectra were obtained for ail other engine 
conditions, but have not yet been analyzed. The HITRAN96 data compilation, soon to be released, 
will contain a separate file of high-temperalure C0 2 lines which may allow a full set of C0 2 
emission indices to be derived from our spectra, but until that analysis is carried out, we cannot be 
certain how accurate that database will prove to be, and whether laboratory characterization of 
unknown lines will also be required. 'Hie C0 2 emission index listed in Table D-2 differs by less than 
1 0 percent from the theoretical value of 3260. The theoretical values were determined from the fuel 
analysis data, assuming that all the carbon in the fuel was converted to C0 2 during combustion. 
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Figure D-2. Schematic of TDL dL^raostic beam path through engine test celh 
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Figure D-4. Reference cell intensity scan (corresponding to Fig. D-3) 
showing reference lines of target gases. 
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Table D-l. Emission Index Values for NO Determined from IR TDL Absorption in the 
1897-cnr' Region, with Total Error Limits and Quantities Used to 
Estimate Them 


Condition 

T3/AItit»dc 

Elt'NO) 

Error 
Limit 
jn Ef 

Standard 
Deviation in 
Optical Depth 

T-Katio EX 
Conversion 
Factor 

t;v lii(NO) 

Conversion 

Factor 

Total Mass 
Flux Factor: 
InletAlodel 

858/7.6 km 

27.38 

5.27 

0.0082 

163.70 

156.35 

1.08 

764/7/, km 

13.28 

2.56 

0.0155 

169,32 

164.71 

1.23 

733/7.6 km 

11.15 

2.16 

0.0219 

178.29 

175 71 

1.2! 

668/7.6 km 

8.9 i 

1.73 

0.0304 

226.21 

209.43 

).!4 

512/9.1 km 

206 

0.47 

01240 

364.26 

249.58 

0.99 

733/9.1 km 

12.55 

2.46 

0.0370 

168.77 

167 79 

1.08 

765/9.1 km 

14.26 

2.75 

0.0153 

161.94 

162.95 

1.08 

808/9.1 km 

18.53 

3.57 

0 0118 

155.77 

156 75 

MO 

67 i/9. 1 km 

9.02 

1.74 

0.0130 

18628 

187.15 

1.22 

727/12.2 km 

10.20 

2.03 

0.0505 

187.87 

186.75 

1.08 

766/12.2 km 

11.85 

2.28 

0.0075 

180.36 

182.43 

) ,09 

665/12.2 km 

7.20 

1.39 

0.0099 

207.91 

203.47 

1.17 

557/12.2 km 

3.85 

0.75 

0.0312 

306.25 

255.78 

1.08 

635/15.2 km 

6 01 

1.22 

0.0647 

257.32 

175.07 

1.15 

671/15.2 km 

8.73 

1.70 

0.0303 

267.98 

217.14 

1.17 

732/15.2 km 

9.8X 

1.90 

0.0061 

238.22 

mm i 

i.n 

764/15.2 km 

10.29 

1.98 

0.0155 

238.31 

254.68 

1.07 

733/SL S 

17,19 

1 

3.33 

0.0219 

145.04 

110.62 

1.17 

660/S LS 

13.3) 

2.57 

0.0193 

169.28 

130.77 

1.16 

651/SLS 

12.59 

2.43 ! 

0.0205 

182.59 

140.46 

1.13 

526/S LS 

3.59 

0.90 

0.1629 

392.55 

277.70 

0.99 

517/SLS 

3.43 

1 26 

0.3)47 

407.16 

307.40 

0.97 

762/51 .S 

18.94 

3.67 

0.0258 

139.86 

11834 

1.19 

861/SLS 

41.78 

8.20 

0.0392 

150.29 

119.17 

i.n 
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Tabic D-2. Emission Index Values Obtained from IR TDL Absorption Spectra 
(Spectral Regions of Lines Used in Wavenumbers, 

Given in Parentheses) 
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APPENDIX E 

AEROSOL SAMPLING MEASUREMENTS (INCLUDING CO, DIAGNOSTIC TEST) 


D. E. Hagen, P. D. Whiteficld, A. R. Hopkins, 1 D. Paladino, M. B. Tiueblood 
Cloud and Aerosol Sciences Laboratory 
University of Missouri-Rolla 
Rolla, MO 65409-0430 USA 

E-1.0 INTRODUCTION 

There is a real need to improve the industry standard for aerosol characterization in jet engine 
emissions. The current "smoke number" technique, although useful, docs not provide the essential 
characteristics: total aerosol and non-volatile aerosol concentration, size distribution, and hydration 
or growth as measured by soluble mass fraction and critical supersaturation spectrum. T hesc 
characteristics are needed in order to effectively monitor the environmental impact of engine 
emissions. They also represent the input parameters required by models designed to describe the 
production and dispersion of engine aercrol cmittants from the combustor, through the exhaust 
nozzle plane to their eventual equilibration with the ambient atmosphere. In this appendix, we 
describe the application of a mobile aerosol monitoring facility, the University of Missouri-Rolla 
Mobile Aerosol Sampling System (UMR MASS), to characterize engine aerosol emissions in 
partial fulfillment of the altitude chamber test, which is the subject of this report. The multi- 
configurational MASS has been employed in both ground and airborne field operations. It has been 
successfully flown on the NCAR Sabreliner research aircraft (Ref. F.-l), the DLR (German 
Aerospace Research Establishment) Falcon research aircraft (Refs. E-2 and E-3), and on the NASA 
DC8 research platform in project SUCCESS. In ground tests, the MASS lias participated in 
numerous jet engine related ground tests at NASA Lewis, Pratt and Whitney. AEDC and MDAE 
(Refs. E-4 and E-5). and has been deployed to the Air Force's Phillips Laboratory to resolve aerosol 
generation problems in a nigh power chemical laser system (Ref. L-.-6). In all cases the measure- 
ments were made on samples taken from a harsh physical and chemical environment: with both 
high and low temperature and pressure, and in the presence of highly reactive gases. 

E-2.0 EXPERIMENTAL APPROACH 

A schematic diagram of the UMR MASS (Mobile Aerosol Sampling System) configured 
for the AEDC Current Advanced Commercial Engine Simulation Emissions Studies is shown in 
Fie. E-l . For this test, the MASS was divided into two work stations. The first, located close to 
the extractive sampling probe selection manifolds, was used in conjunction with the extractive 
sampling system to acquire total condensation nuclei (CN) concentrations for all aerosols in real- 
time and to fill and pressurize sample tanks. The exhaust sample was delivered to the MASS from 
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individual and ganged probe sources. An isobaric cooler for supersaturation pulse simulation, a 
laser aerosol spectrometer for real-time large diameter aerosol measurements (diameters from l .0 
to 30 pm), and a needle-to-grid electric precipitator for aerosol collection on electron microscope 
grids were also located at the first work station. The second work station consisted of two fcAC's 
(Electrostatic Aerosol Classifier) configured in series as shown in Fig. E-l . With this system the 
tank samples; taken at the first work station were analyzed for size distribution, and growth and/ 
or hydration (critical supersaturation) information for aerosols in the submicron range 10 nm to 
0.5 pm. 

The MASS methodology is described in the following sections. 

E-2.1 TOTAL CN CONCENTRATIONS 

For total CN the incoming sample stream is routed straight to aCN counter (CNC). Aerosol 
concentrations are measured continuously with typically a 3-sec counting period. In situations 
where the sample pressure drops below 0.5 atm, it is found that particle counting efficiency falls off 
and must be accounted for. The efficiencies for the counters used in this work were measured during 
a recent study devoted to calibration and operation of CNC's. 

E-2.2 SIZE DISTRIBUTION MEASUREMENT 

Aerosol size distribution measurements are made using a single EAC (Ref E-7) where 
sample aerosol laden air is drawn from the source and ducted to the MASS (see Fig. B-l ), In cases 
where the sample is available only for a short time (i.e., seconds), an intermediate storage tank is 
used to store the sample for subsequent processing with the MASS. For size distribution and hydra- 
tion characterization, the air is directed along path (1) and given an electrical charge in the bipolar 
charger. The aerosol then passes through HAC1 and a selected size aerosol (monodispersed) 
emerges. Chian sheath air is provided by the pump and absolute filter. The monodispersed aerosol 
is directed around EAC2 via path (2) and on into the CN counter. EAC1 is stepped through a 
selection of voltages (aerosol particle sizes) and a particle concentration is taken at each particle 
size with the CN counter, and hence a size distribution is generated. The aerosol can be directed 
around both EAC’ s via path (3) and sent directly to the CN counter to get a total CN count. The 
count versus voltage data set must be numerically analyzed in order to generate the actual size 
distribution (Ref ESr8). Hero we note that an interesting feature of the EAC is (hat it not only takes 
a polydisperse aerosol and separates out a monodisperse fraction, but it also transfers this 
monodisperse fraction to the clean sheath air stream. The EAC can likewise be used to remove 
aerosol from the original sample ai: stream. This methodology applies for both total aerosol 
sampling and for non-volatile aerosol sampling. In the latter case an in-line preheater set at 175’C 
provides the means to separate the non-volatile component of the total aerosol flow. 
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E-2.3 AEROSOL HYDRATION EXPERIMENT 

The UMR MASS can access aerosol hydration properties through the haze chamber or 
Laktionov technique (Refs. E-8, 9, and 10). Sample air is drawn into the MASS and size selected 
in EAC1 as above. This preselected monodispersed aerosol passes via path (4) to the saturator 
integrated into EAC2. Here it is brought to 100-percent relative humidity in the saturator, and the 
aerosol particles deliquesce to their equilibrium size at 100-percent relative humidity. Even though 
the particles entering the saturator are monodispersed in size, a spectrum of particle sizes will exit 
the saturator because of differences in particle chemistry and soluble mass fraction. To investigate 
the deliquesced particle size spectrum, the voltage on EAC1 is held fixed. The voltage (size setting) 
on E AC2 is then stepped through various values. The aerosol is passed onto the CN counter, and a 
count is taken for each voltage setting of EAC2. These counts can then be translated into a critical 
supersatu ration spectmm for the particles of the size given by EAC1 . 

Under conditions where particles are sampled at or near the exit plane of an exhaust nozzle, 
the particle emission index can be determined by measurement of the number of particles emanat- 
ing from the engine per cm 1 of gas in the exhaust flow, and knowledge of fuel/nir ratio at the sample 
location. The latter quantity can be simply related to the mass of fuel (in kilograms) for each cm’ of 
gas passing through the engine. 'Hie ratio of the number of particles per cm’ to the kilograms of fuel 
per cm' yields the number-based emission index. A mass-based El has also been estimated (Table 
E-2). In this case a soot volume fraction (SVF) was extracted from appropriate size distribution 
data, assuming all particles to be spherical and of known diameters. This SVF was then converted 
into mass units by employing an average particle density of 2 gm/cnr 5 . 

E-3.0 PARTICLE CHARACTERIZATION RESULTS 

In this project the following aerosol properties were characterized as a function of altitude. 
and T3 and P3: total concentration, aerosol, number- and mass- based emission indices, soot volume 
fraction, total aerosol and non-volatile aerosol size distributions and hydraiion (soluble mass frac- 
tion). Table F.-l is a compilation of the aerosol concentration [CN] and number-based emission in- 
dices El(CN). Information relating to data acquisition start times, test facility data acquisition 
steady-state data point (SSDP) number, sample source defined by probes sampled, T3, P3, percent 
CO-, and altitude are also provided. Figures E-2 through F.-4 are correlations made between these 
data and other test parameters, i.e., Fig. E-2 is a plot of the aerosol number-based emission indices 
as a function of thrust (a*, represented by CO, concentration) at simulated sea-lcvcl-static (SLS); 
Fig. E-3 is a plot of aerosol number-based emission indices as a function of altitude; and Fig. E-4 
is a plot of aerosol number-based emission indices as a function of T3 for 9.1 km. These data show 
no strong correlation between aerosol emission index and thrust, altitude, and combustor inlet tem- 
perature. However, there appears to be a correlation between combustion efficiency and aerosol 
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emission index observed in the ground idle data points, where elevated emission indices for aero- 
sols, CO, and unburned hydrocarbons are recorded. 

Size distributions for total aerosol and non-volatile aerosol samples were recorded for a 
selected subset of conditions and sample sources. Table E-2 lists the time, test facility SSDP 
numt>cr, sample source, and UMR data point (run number) for all of the size distributions measured. 
Piots of these distributions arc given in Figs. E-5 through E-14. In these figures each size 
distribution is identified by UMR data run number, time, test facility SSDP, and altitude. 

Hydration data analyzed as the fraction of a given aerosol sample population as a function of 
soluble mass fraction were acquired for aerosol dry sizes 20. 50, and 80 nm and continue to be 
deconvolved. The general observation from these data is that very little soluble mass is found on 
these aerosols with very few aerosols having a soluble mass fraction > 0.1 . The dry sizes of 20, 50, 
and 80 nm were chosen because these sizes tend to encompass the peak in the size distribution 
typically observed for aircraft engine soot emissions. 

E-4.0 CONCLUSIONS 

The purpose of this report is to provide the scientific community with the data from this 
altitude chamber test as soon as possible. As a result at this time only the following preliminary 
conclusion 4 : can be drawn from these data. This is the first database of its kind and provides 
particulate characterization (i.e., total concentration, aerosol, number- and mass-based emission 
indices, aerosol (soot) volume fraction, total aerosol, and non-volatile aerosol size distributions and 
hydration (soluble mass fraction)) for an advanced current commercial engine simulation operating 
at: 0.8 Mach cruise conditions at altitudes from simulated SLS to 15.2 km. These data characterize 
the aerosol generated in the engine and that continues to exist at a distance <10 cm downstream of 
the exit plane of the engine exhaust nozzle. 

There is no strong correlation between emission index and thrust (measured as CG 2 
concentration), altitude, or combustor inlet temperature (T3). However, there is a correlation 
between combustion efficiency and aerosol emission index. In the case at SLS where CO and 
unbumed hydrocarbon emission indices are elevated, a similar elevation is observed in the aerosol 
emission index. The mean aerosol number-based emission index calculated from all measurements 
L<; (2.2 ± 0.7) x 10 l \ the associated mean mass-based emission index is 0.012 ± 0.001 (where the 
quoted uncertainty in both cases is the standard deviation of the mean). 

Typical sizii distributions arc log-normal and peak in the vicinity of 20 to 40 nm diameter. 
The volatile component is small and in most cases can be considered negligible. 


127 


AEOC-TR-96-3 


When exposed to a supersaturation pulse, the aerosol concentrations in sample flows were 
enhanced by at least a factor of four. 

E-5.0 C0 2 DIAGNOSTIC TEST 

The McDonnell Douglas CO, detector was designed to be an efficient, compact, and portable 
means of detecting CO, to be used in conjunction with the Mobile Aerosol Sampling System 
(MASS) in order to obtain emission indices during ground testing of aircraft engines. 

Air from the engine exhaust flowed through two absoiption cells before being exhausted to 
the mechanical purnp (Fig. E-l 5). The flow was regulated with a flowmeter located downstream of 
the cells. The flow meter was set to I liter/min. The pressure in the cells was continually monitored 
using a capacitance manometer. 


The pressure in the absorption cells was maintained near 1 atm by regulating the flow through 
the detector with a valve located upstream of the CO, detector. The air from the engine exhaust was 
preconditioned before entering the detector by flowing the air over Dryerile* to remove water vapor 
and by flowing it through a filter to remove particles. 

The absorption cells are commercially obtained units which detect absorption of the v ; , 
band of C0 2 near 4.25 pm. These cells were connected in parallel to each other and to the 
capacitance manometer. 

The absorption cells each consisted of a modulated infrared light source, an infrared filter 
centered at 4.2.5 pm, an absorption cell, and a detector. The path lengths of the two cells were 
fabricated to maximize sensitivity for Cell 1 of 0 to 2 percent CO, and 0 to 3000 ppm for Cell 2. 

E-S.l CALIBRATION VERSUS KNOWN CO, CONCENTRATIONS 

The absorption cells are calibrated at the factory to read in percent CO, for a gas at standard 
temperature and pressure. The absorption of the light is determined by Reel ’s Law, log(l/b) = ~cL, 
where L is the unatieriuatcd light intensity, l is the attenuated intensity, L is the length of the cell, 
and e is the absorption of CO,. The: absorption. £, is proportional to the density of CO, in the cell; 
therefore, the cells are calibrated at the factory with gases of known percentages of CO, in the range 
of interest. Calibration curves are then generated over the range of interest and fit with a 
mnltiparametcr fitting equation. These curves are then programmed into a microprocessor so that 
the microprocessor outputs a voltage that is linearly proportional to the percent of CO, for a given 
detector attenuation (I/L). This process is accurate only over a limited range of CO, partial pressures. 
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Wc have increased this range by employing two absorption cells of differing ranges. In this 
manner we can monitor C0 2 percentages from 10 percent down to below 300 ppm. We checked the 
factory calibrations in our laboratory using gases with known concentrations of C0 2 . These 
included calibration gases of 2 percent, 0.2 percent, 1000 ppm and 300 ppm. 

E-5.2 MEASUREMENT OF THE ENGINE EXHAUST 

The purpose of employing our detector in these tests was to validate our detector by 
comparison with standard techniques. The C0 2 in engine exhaust was measured by employing the 
procedure described above. C0 2 concentrations were derived using the readings from the 
absorption cells. Absolute densities of C0 2 were determined using the CO } readings and the 
pressure and temperature of the absorption cells and the ideal gas tew (PV » nRT). 

The system as currently constructed is limited to making measurements at cell pressures at or 
above ambient pressure. This is the normal condition encountered under ground test conditions. 
Therefore, we limited our comparisons with other measurements taken by other groups during the 
engine testing to those which were taken at above ambient pressures. These comparisons showed 
our measurements to be always within 5 percent of the AEDC C0 2 measurements using an NDIR 
analyzer (Appendix B). The AEDC CC) 2 gas sampling measurement results in turn compared 
favorably with the fuel-to-air ratios determined from the test facility engine performance data. 

The system can be enhanced to make measurements at sub-ambient pressures by employing 
a two-chamber vacuum pump so that the blank-off pressure of the system is much less than the 
pressure in the engine exhaust. 
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Figure E-l. UMR Mobile Aerosol Sampling System. 
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a. T3 = 733 K, Altitude = 9.1 km. 
Time - 09:34:00 


b. T3 = 733 K, Attitude - 9.1 km, 
Time = 10:27:20 



c.. T3 = 733 K, Altitude - 9.1 km, 
Time = 10:29:00 



d. T3 * 733 K, Altitude = 9.1 km, 
Time = 10:42:08 



e. T3 = 733 K, Altitude ~ 9.1 km, 
Time a 10:47:00 


Figure K-6. Total aerosol size distribution, Set 2. 
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a. T3 = 727 K, Altitude ~ 12.2 km, 
Time* 12:07:00 


b. T3 = 727 K, Altitude - 12.2 km, 
Time = 12:28:00 




c. T3 * 7(66 K, Altitude * 12.2 km, 
Time = 12:52:00 


d. T3 = 766 K, Altitude = 12.2 km, 
Time * 12:56:00 




e. T3 * 766 K, Altitude = 12.2 ktn, 
Time * 13:00:00 


f. T3 ~ 66S K, Altitude * 12.2 km, 
Time = 13:20:00 


Figure E-8. Total aerosol size distribution, Set 4. 
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Digmrtitr, cm OUraeUr.cm 


a* T3 = 733 K, Altitude = 3,1 km(SLS), *>. T3 = 733 K, Altitude ~ 3.1 km{SLS), 
Time - 15:05:00 Time « 15:38:00 



2 9 « 5 t 7 8 2 10 * 10 * 

10** 10* OltWtt*f.om 

Dl»m*t*r. era 


c. T3 = 733 K, Altitude » 3.1 km(SLS), d. T3 = 669 K, Altitude » 3.1 km(SLS), 
Time ~ 16:07:00 Time ~ 16:24:00 



e. T3 = 651 K, Altitude =: 3.1 km(SLS). f. T3 = 526 K, Altitude « 3.1 km(SLS), 
Time - 16:36:00 Time * 16:45:00 

Figure E-10. Total aerosol size distribution, Set 6. 
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1 1 4 $ 6 7 « , 2 

i<f* to 4 

DUin#t*?, cm 


a. T3 and Altitude unknown, 
Time = 07:57:00 



Climate', cm 


c. T3 a 733 K, Altitude - 7.6 km, 
lime = 08:39:50 



2 J 4 5 « r 8 


JO* 10' 1 

M tender, cm 

b, T3 - 764 K, Altitude * 7.6 km, 
Time = 08:22:34 



2 $ 4«47» 

1& 4 10 * 


d. T3 and Altitude unknown, 
Time = 08:59:30 



0j»rr*4*f, cm 


e. T3 = 733 K, Altitude = 9.1 km, 
Time = 10:29:00 


Figure E-12. Non-volatile aerosol size distribution, Set 5. 
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a. T3 = 733 K, Altitude = 3.1 km(SLS), 
Time = 15:05:00 



b. T3 = 669 K, Altitude = 3.1 km(SLS), 
Time = 16:24:00 



c. T3 = 651 K, Altitude = 3.1 km(SLS), 

Time = 16:36:00 

Figure E-14. Non-volatile aerosol size distribution, Set 3. 
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Tnble E-l- CN Concentration and Number-Based El as a Function of Test Parameters 
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T3, K 
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[CN] 

%CO a 

El (CN) 

AH, km 
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6:26 14 
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60000 0 
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7.0 

31410 

8:43:30 
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9.426411 
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9.1 
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5.35E+11 
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11.872 
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9:31 :00 

37.38,40-48 
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34320 
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37.38,40-49 
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3.37 
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0.1 
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11.872 
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91 
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P5 
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P4 
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11.072 
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40560 
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3.03 
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9.1 

40620 

1 1:17:00 
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3.73 
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7.80E+12 

9.1 
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Table E-l. Continued 


Tkm t 

sec 

Time 

Facility 

8$DP 

Sample Source 

T3.K 

P3, atm 

EC fel 

<KCO t 

Cl (CN) 

AK t km 

41370 

11:29:30 

54,55 

160 
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14.576 

127576.9 

4.11 

5.79E+12 

91 

41460 

11:31:00 

54,55 

360 
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14.576 

120059.0 

401 

5.61E+12 

9.1 

42180 
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57,56 
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8.779 
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91 
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2.70 
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2.77 
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3.53 
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3.40 
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12.2 
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3.54 
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3.53 
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3.40 
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12 2 
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8.300 
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3.54 
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727 
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12.2 
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102085.4 
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12.2 
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12.2 
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706 
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12.2 

46440 

12:54:00 

66,67,68 

180 

786 

9.099 

1939612 

3.70 

1.57E+13 

12.2 

46500 

12:55:00 

66,87,66 

360,180 

786 

9.099 

170380.9 

3.70 

1.23S+13 

12.2 

46620 

12:57:00 

68.67.68 

360 

766 

9.099 

146296.6 

3.70 

1.03£*13 

12.2 

46700 

12:58:20 

68,67,68 

180 

766 

9.099 

199163.9 

3.70 

1.39E+13 

12.2 

46760 

12:59:20 

66,87,88 

360,180 

766 

9.090 

163662.2 

3.70 

1.14E+13 

12.2 

47960 

13:19:20 

71.72 

360,180 

065 

0,331 

50550.8 

2.60 

0.34E+12 

12.2 

48060 

13:21:00 

71,72 

360 

665 

6.331 

21069.6 

2.86 

2.34E*12 

12.2 

48120 

13:22:00 

71,72 

180 

665 

0.331 

20870.6 

2.88 

2.22E+12 

12.2 

46160 

13:23:00 

71.72 

360,180 

665 

0.331 

24348.0 

2.98 

2.59E+12 

12.2 

50330 

13:58:50 

79 

180,360 

071 

4.011 

8900.0 

3.16 

1.52E+12 

15.2 

51420 

14:17:00 

63,84 

180,360 

764 

5.045 

340000.0 

3.98 

3.79E+13 

15.3 

53610 

14:53:30 

86-94 

P10 

733 

10.741 

223584.4 

3,10 

135E+13 

3.KSLS) 

53640 

14:54:00 

86-84 

P11 

733 

18.741 

212186.0 

3.29 

1.23E*13 

3.1 (SLS) 

53700 

14:55:00 

86-94 

P12 

733 

18.741 

198167.0 

3.25 

1.15E+13 

3.1 (SLS) 

53760 

14:56:00 

86-94 

P13 

733 

10.741 

184831.7 

2.43 

1.44E*-t3 

3.1 (SLS) 

53880 

14:58:00 

86-94 

P14 

733 

10.741 

156344.3 

1.39 

2.10E+13 

3.1 (SLS) 

83950 

14:59:10 

88*94 

P33 

733 

10.741 

97618.1 

1.39 

130EM3 

3.1 (SLS) 

54020 

15:00:20 

56-94 

P32 

733 

10.741 

211030.0 

2.43 

1.07E+13 

3.1 (8LS) 

54060 

15:01:00 

86-94 

P31 

733 

10.741 

230730.0 

3.25 

1.37E+13 

3.1 (SLS) 
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Time, 

MS 

54120 

54130 

54250 

54375 

54420 

54460 

57280 

57320 

57300 

57420 

57550 

57510 

57675 

57750 

57780 

57840 

57000 

67960 

59100 

59850 

59880 

59940 

60360 

60420 

60980 

61020 

61620 

61740 

61800 

32260 

62400 

62460 


Table E-l. Concluded 


Time 

Facility 




SSDP 

$ampt« Source 

T3.K 

15:02:00 

56 '4 

P30 

733 

15:03:00 

86-94 

P28 

733 

15:04:10 

86*94 

P2B 

733 

15:06:15 

66*94 

180 

733 

15:07:00 

86-94 

360 

733 

15:06:00 

86-94 

180,360 

733 

15:54:40 

86-94 

P19 

733 

15:55:20 

98-94 

P20 

733 

15:56:00 

86-94 

P21 

733 

15:57:00 

86-94 

P22 

733 

15:59:10 

86-94 

P23 

733 

16:00:10 

66-94 

P24 

733 

16:01:15 

88-94 

P6 

733 

< S8:02:30 

86*94 

P5 

733 

16:03:00 

86-94 

P4 

733 

16:04:00 

86-94 

P3 

733 

18:005 

86-94 

P2 

733 

16:06:00 

86-94 

PI 

733 

16:25:00 

86,96 

380 

669 

16:37:30 

100.10? 

360 

651 

16:38:00 

100,10! 

180 

651 

16:39:00 

100.101 

360,180 

651 

16:46:00 

103,104 

380 

526 

16:47:00 

103,104 

160 360,180 

528 

16:56:00 

105,106 

3,0 

517 

16:5700 

105,106 

160 

517 

17:07:00 

106,109 

360,180 

762 

17:09:00 

108,109 

360 

762 

17:10:00 

108.109 

180 

762 

17:18:00 

lit, 112.113 

360,160 

861 

17:20:00 

111,112,113 

360 

861 

17:21:00 

111.112,113 

160 

861 


P3, atm 

fCNl 

%co 2 

El (CM) 

Aft, km 

16 741 

238004.6 

3.29 

1.37E+13 

3.1 (SLS) 

18.741 

220117.0 

3.16 

1.33E+13 

3.1 (SIS) 

18.741 

188277.1 

3.09 

1.18E+13 

3.1{$LS) 

18.741 

235153. S 

3.09 

1.45E+13 

3.1 {SIS) 

13.741 

2170297 

312 

1.33E+13 

3.1 (SLS) 

16741 

2103156 

3.13 

1.27E+13 

3.1 (SLS) 

18 741 

236441.9 

3.00 

1.92E+13 

3.1 (SLS) 

18741 

263850.1 

3.14 

1.60E+13 

3 USLS) 

18741 

256059.3 

328 

1.49E+13 

3.1 (SLS) 

18741 

244430.3 

3.27 

1.42E+13 

3.1 (SLS) 

18741 

210031.0 

2.17 

1.84E+13 

3.1 (SLS) 

18741 

142027.9 

0.82 

3.30E+13 

3.1 (SLS) 

16741 

165384.1 

1.03 

3.06E+13 

3.1 (SLS) 

18741 

164247.0 

2.12 

1.49E+13 

3.1 (SLS) 

16741 

196426.0 

3.07 

1.23E+13 

3.1(SLS) 

18741 

212583 2 

3.24 

1.25E+13 

3.1 (SLS) 

18.741 

253221.3 

3 17 

1.56E+13 

3.1 (SLS) 

16741 

190684.3 

3.37 

1.10E+13 

3.1 (SLS) 

13.330 

391127.1 

2.54 

2 99E-M3 

3.1{SLS) 

12041 

279557.7 

2.43 

2.28E+13 

3.1 (SLS) 

12.041 

229123.8 

2.41 

1.89E+13 

3.1 (SLS) 

12,041 

189376 1 

2.44 

1.49E+13 

3.1 (SLS) 

5.158 

34077.3 

1.47 

439E+12 

3.1 (SLS) 

5.158 

15605.2 

1.48 

2.43E+12 

31 (SLS) 

4.768 

507454.3 

1.42 

8.20E+13 

3.1 (SLS) 

4.766 

437381 1 

1.41 

7.20E+13 

3.1 (SLS) 

21.473 

328168.1 

3.33 

1.81E+13 

3.1 (SLS) 

21.473 

417273.6 

3.32 

2.31E+13 

3.1 (SLS) 

21.473 

273896 7 

3.20 

1.56E+13 

3.1 (SLS) 

30.458 

2504 75.3 

4.29 

1.08E*13 

3.1 (SLS) 

30.458 

378987 5 

4.28 

1.046+13 

31 (SLS) 

30.458 

315494 3 

424 

1.38E+13 

3 1 (SLS) 
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APPENDIX F 

OH-UV LASER ABSORPTION SYSTEM 

C. W. Brasier 
LA. Drakes 
R. P. Howard 

F-I.O INTRODUCTION 

A laser absorption system was installed in the test cell during the NASA test sequence as part 
of an ongoing effort to increase the diagnostic capability of ABDC. This effort represented the first 
turbine test cell entry of the AEDC UV laser measurement system. Consequently, the majority of 
the effort was aimed at understanding the dynamics of the test cell interactions with the laser 
system, even though a high priority was placed on (he acquisition of data on OH concentration lev- 
els during this engine test period. Although OH absorption was not directly observed in the exhaust 
flow, the data are useful in that upper bounds to the OH concentration may be discerned from the 
detection threshold. 


F-2.0 METHOD 

The procedure for conducting and interpreting an absorption measurement is well 
established. For this application, a primary assumption was that the flow field was spatially uniform 
across the entire exit plane. In that case, the number density, N„ of the OH in a given rotational state 
is determined by the absorption out of that state, denoted as Mo an upper state, denoted as 2 , via, 


= 



( 1 - 1 ) 


where g, and g, represent the state degeneracies, A 2J is the Einstein coefficient for the transition, v„ 
is the frequency of the transition, and k r is the absorption coefficient given by 


t, = tfa('f) - ifo(i) 


(F-2) 


in this expression* /„ refers 10 the measured incident laser intensity, / refers to the measured 
laser intensity after passing through the media. L refers to the physical path length of the media, and 
t v is the spectral transmission. 


ISO 
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For this application, no absorption profiles were observed, so an assessment of the threshold 
detection was made. First, the absorption data from the individual laser frequency scans were aver- 
• aged together and the nominal peak-to-peak noise was noted. To be conservative, the peak-to-peak 

noise variation was considered to be the largest possible signal due to OH absorption which could 
be present in the data. Id other words. 


k v<i L<b pp (k v L) 


(F-3) 


The peak height of the absorption line can be expressed as 


i - s F 2 e * - 1 " er f ( a } } 

v o ~ 4 n~ X D 


(F-4) 


where S is referred to as the line strength, which is simply the frequency integral of Eq. (F-2). The 
lint! strength can be directly converted into the state specific number density of the lower, absorbing 
state, N„ and the Einstein coefficient of the transitions (i.e., A }l ), through Eq. (F-l). For this work, 
the Einstein coefficients were taken from first-principle calculations. In Eq. (F-4). a is the broaden- 
ing parameter, which was taken as 0.2 for this work. Note that the results are extremely insensitive 
to the value of a chosen. The remaining parameter is the Doppler width of the spectral line, y c , 
which is dependent on the gas sialic temperature. The appropriate reported temperature measure- 
ments from the probe data were used to determine Furthermore, the characteristic path length, 
L. was taken from the half-height points of the pressure profiles. Thus, by combining these factors 
with Eq. (F-3), it was possible to define an upper limit to the state specific number density. 

In order to determine the total OH density from the state specific upper limit, the partition 
function was evaluated under the guiding assumption that the rotational and translational modes 
were in at least near equilibrium with each other, implying that the above probe temperatures could 
be used. The total OH(X) density then follows as 


N| Q(T) 

OH exp (-£,/* 7')’ 


(F-5) 


where Q(T) is the temperature-dependent partition function. 
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F- 3.0 APPARATUS 

Experimentally, the measurement followed the standard procedure for an amplitude- • 

modulated absorption measurement, with a few modifications specific to this application. Shown 
in Fig. F-l ss a schematic of the experimental setup used for the OH absorption measurements. 

The system consisted of a continuous wave Coherent 699-21 scanning ring dye laser pumped by 
a Coherent INNOVA-400 argon-ion laser. Using Kiton Red laser dye, an intra-cavity Lithium 
Niobatc frequency doubling crystal, and a pump laser power of 8 W at a wavelength of 514.5 nm, 
the ring laser output was nominally 3 mW single frequency at a wavelength of 306 nm. The 
frequency -doubled dye laser wavelength could be tuned over a wavelength range of approximate- 
ly 304 nm to 327 nm with a frequency iinewidth on the order of 500 kHz, which corresponds to 
a wavelength resolution of roughly 1 .5 x 10-" nm at 306 nm. The dye laser output frequency could 
be scanned ± 30 GHz about the selected absorption line of the OH molecule. A reference signal 
from the laser controller representing the laser frequency scan range was recorded bv the data 
acquisition system. 

A Burleigh wavemetor, model WA-10, and an OH discharge tube were included in the 
setup to provide laser diagnostic instrumentation. The wavemeter provides a convenient method 
of monitoring the fundamental laser wavelength, centered about 6 1 2 nm, during initial alignment 
and laser wavelength selection with a wavelength resolution of 0.003 nm The frequency-doubled 
output of the dye laser was directed onto a retractable mirror which, when placed in the beam 
path, directs the beam to the OH discharge tube for additional verification of the proper laser 
wavelength. Once the OH absorption feature from the OH discharge tube was verified, the 
retractable mirror was removed, and the beam was then passed through a chopper wheel 
operating at 1.4 kHz and focused onto a 30-m, 200-jim core diameter (JV grade multimode fiber 
optic cable for delivery to the test cell. 

The beam arrangement in the test cell is shown schematically in Fig. F 2. At the fiber optic 
output the beam was split. Approximately 5 percent of the beam was directed onto a Hamamatsu 
I P28 photomultiplier tube with a UV-cnhanccd photocathode for recording the laser power exiting 
the fiber optic cable. The photomultiplier output was passed through a pre-amp and a Stanford 
Research, model SR-250 lock-in amplifier referenced to the 1.4 kHz chopper frequency and 
recorded in the data acquisition system. This reference signal was remotely monitored from the 
laser room so that minor adjustments to the kixer beam alignment as it entered the fiber optic cable 
could be made in order to optimize the throughput. 

After passing through the beam splitter in the lest cell, the majority of the beam was directed 
across the turbine engine exit plane and onto n Hamamatsu 1P28 photomultiplier tube with UV- 
enhariced photocathode and pre-amp to provide the OH absorption measurement. The pre-amp out- 
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put was directed to a second Stanford Research model SR-250 lock-in amplifier referenced to the 
1.4-kHz chopper frequency and recorded by the data acquisition system. 

The data acquisition system consisted of a Data Translation 16-channe! 12-bit A/D board 
inside a microcomputer. The laser frequency scan width, the fiber optic reference power, and the 
absorption signal were all recorded by this unit. During the OH measurement the laser frequency 
scanned a frequency interval of ± 30 GHz in 2.5 sec per scan. The data from the reference photo- 
multiplier tube, the absorption photomultiplier tube, and the laser scan width were each sampled at 
l .0 kHz during the data acquisition period. 

F-4.0 RESULTS 

One result was glaringly evident from this measurement effort: the OH density at the exit 
plane was below the sensitivity threshold of the absorption system. System cross-checks, such as 
verifying the wavelength setting with the OH discharge lamp, during these measurements indicated 
that the entire system was functioning properly. Yet, no discernible trace of OH absorption was 
detected even though the strong Q, lines were monitored. This was true regardless of the altitude 
and engine power settings throughout the test period. 

An effort was made to estimate the upper bound of the OH number density for the higher 
power settings of the 9.1 -tin altitude using the algorithm described above. The calculational results 
are shown in Table F-l. From the variety of lines measured, it is clear that the upper bound to the 
maximum OH number density is nominally 
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Table F-l. Intimates of Maximum OH Number Density at the 9.1-km Condition 


T3K 

Line Measured 

Maximum OH, cm” 1 

671 

Qi( 5) 

1E12 

733 

0,(2) 

IEI2 


0,(3) 

2E12 


Q;( 4) 

2E12 


0.(5) 

2E12 


0.(6) 

1E1? 

765 

Q;(4) 

3E11 

803 

Q,(4) 

5E1 1 


Thorne A P . S}>at> '>ph\ju'f. Chapman ant! Hit!!. 1974. 
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